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Abstract 



The potential between a static quark and antiquark in pure SU(3) Yang- Mills theory 
is evaluated non-perturbatively through computations on the lattice in the region 
from short to intermediate distances (O.OSfm < r < 0.8fm). Renormalized dimen- 
sionless quantity are extrapolated to the continuum limit, confirming the theoretical 
expectation that the leading lattice artifacts are quadratic in the lattice spacing. 
In the high energy regime the results are compared with the parameter-free predic- 
tion of perturbation theory obtained by solving the Renormalization Group equation 
at two and three loops. The choice of the renormalization scheme to define a running 
coupling turns out to be important for the accuracy of the perturbative prediction: 
by obtaining the running coupling through the force, perturbation theory is applica- 
ble up to q; ~ 0.3, while from the static potential only up to a ~ 0.15. In the region 
where perturbation is supposed to be reliable, no large unexpected non-perturbative 
term is observed: on the contrary, one finds a good agreement between perturbation 
theory and our non-perturbative data. 

For large quark-antiquark separations our results are compared with the predictions 
of a bosonic effective string theory, finding a surprising good agreement already for 
distances > 0.5fm. 

In the second part of this work, universality and scaling behavior of different formu- 
lations of Yang-Mills theory on the lattice are discTissed. 

In particular, the Iwasaki and DBW2 action are investigated, which were obtained 
by following renormalization group (RG) arguments. The length scale Tq ~ 0.5fm 
is evaluated at several lattice spacings and the scaling of the critical deconfinement 
temperature T^Vq is analyzed and confronted with the results obtained with the usual 
Wilson plaquette action. Since they agree in the continuum limit, the universality 
is confirmed. We remark that the quantity to use to set the scale has to be chosen 
with care in order to avoid large systematic uncertainties and Tq turns out to be 
appropriate. For the critical temperature the data obtained with RG actions show 
reduced lattice artifacts, above all with the Iwasaki action. 

Finally the mass of the glueballs O^"*" and 2++ is evaluated by considering the quan- 
tities mo++ro and m2++ro; however for those observables no clear conclusion about 
the scaling behavior can be drawn. Particular attention is dedicated to the viola- 
tion of physical positivity which occur in these actions and the consequences in the 
extraction of physical quantities from Euclidean correlation functions. 
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Introduction 



The strong, weak and electromagnetic interactions are up to now described success- 
fully within the Standard Model of elementary particles. This theory is based on a 
local gauge principle [1] , with the gauge group 

Gioc = SU(3)c X SU(2)^ X U(l)^, (1) 

and it is essentially determined once the matter fields and their transformation laws 
under Gioc are specified. 

The degrees of freedom are respectively the color for SU(3), weak isospin for SU(2) 

and weak hyperchargc for U(l). 

In particular, quantum chromo dynamics (QCD) associated to the color group SU(3) 
is the currently accepted framework to describe strong interactions. The matter 
components are the quarks, which are described by spinor fields carrying three color 
indices and appearing in six flavors (up, down, charm, strange, top and bottom). 
The gluons arc the vector bosons that mediate the interactions and, due to the fact 
that the gauge group SU(3) is non-Abelian, self-interactions exist, contrary to the 
situation in electromagnetism. 

The formulation of QCD was the outcome of several years of investigations and dis- 
cussions on experimental results. 

The goal was to formulate a theory which is able to reproduce the two main features 
of strong interactions, asymptotic freedom and confinement, which are related to 
their properties in the high and low energy regimes, respectively. 
The first step was the observation that the hadrons, which are the particles that 
participate in strong interactions, fall into multiplets that should reflect underlying 
internal symmetries. To express this fact in a concrete way, it was formulated the 
hypothesis that hadrons are composed of more elementary constituents with certain 
basic symmetries, which were called quarks [2,3]. 

Asymptotic freedom was observed in deep inelastic scattering: at high energies 
the form factors in the cross sections loose their dependence on certain mass pa- 
rameters; this scale invariance (also called scaling) was interpreted to mean that the 
constituents act as if they were free particles at extremely high energies. The most in- 
tuitive explanation of the scaling relation came from Feynman's parton model [4-6], 
where the proton was assumed to consist of point-like constituents. Later the par- 
tons were identified with the quarks on the basis of their quantum numbers. 
Finally in 1973 Fritzsch, Gell-Mann and Leutwyler [7] proposed QCD as the non- 
Abelian gauge theory built from the SU(3) gauge group associated with the color 
symmetry. 
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The most direct experimental evidence that quarks must have an additional quan- 
tum number, the color, came from the measurements of the total cross section for 
the annihilation of electron-positron pairs in colliding beam experiments (SLAC, 
PETRA). 

In the same year, it was shown [8-10] that non-Abelian gauge theory exhibits asymp- 
totic freedom. 

The appropriate framework to discuss this property is the renormalization group 
equation. 

The insertion of quantum corrections leads to divergences in the evaluation of phys- 
ical observables; these can be eliminated via renormalizing the theory by adding 
counterterms in the Lagrangian in a proper way. The consequence of the renor- 
malization of a quantum field theory is that the coupling constant that appears in 
the Lagrangian becomes a running coupling. The so-called /3-function quantifies the 
dependence of the coupling on the energy scale. Through the perturbative study of 
the /3-function it was shown that at very high energies (or equivalently at small dis- 
tances) the coupling becomes small, so that the quarks behave as free particles. In 
this regime one then expects that perturbative methods furnish reliable predictions 
for physical observables. 

On the other hand, quarks have never been detected in isolation, but only as 
constituents of hadrons. Prom the theoretical point of view, this should correspond 
to the fact that all physical states are singlets with respect to the color group. In 
order to check whether this feature is contained in QCD, one can not apply pertur- 
bative methods, since the coupling is expected to be large at scales corresponding to 
the size of hadrons. It is however widely believed that this behavior is a consequence 
of quantum chromodynamics, although up to now no proof exists. 
One of the strongest evidence for this to be true comes from lattice field theory. 
In 1974 Wilson proposed [11] a formulation of a gauge field theory on a discretized 
Euclidean space-time. In this framework, by using a strong coupling expansion it 
was possible to demonstrate that at sufficiently strong couplings, pure SU(3) lattice 
gauge theory exhibits color confinement: in the limit of large distances, the energy 
needed to separate a pair of quark- antiquark grows proportionally to the separation. 
Since then lattice field theory has been widely studied and became a powerful tool 
to investigate properties of strong interactions. 

The natural framework to quantize the lattice theory is the path integral formalism; 
in this formulation the system takes the form of a classical four-dimensional statis- 
tical model. The matter fields are treated as classical stochastic variables assigned 
to the points of the lattice, while the gauge fields are associated to the links. In 
this analogy, the Euclidean action of the quantum field theory corresponds to the 
classical Hamiltonian of the statistical system and the mass of the lightest particle 
corresponds to the inverse of the correlation length. 

The lattice furnishes a regularization of the theory by providing an ultraviolet cutoff 
proportional to the inverse lattice spacing and this is actually the only known non- 
perturbative regularization of quantum field theory. 

The crucial question is whether the theory formulated on the lattice is well defined 
in the limit when the lattice spacing a is sent to zero. The physical continuum limit 
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corresponds to holding the values of physical quantities fixed while letting a ^ 0; 
the renormalization group equation describes how the parameters behave by chang- 
ing the scale of the theory (in this case the lattice spacing). It turns out that the 
continuum limit is realized when the bare gauge coupling qq is sent to zero. In the 
language of statistical mechanics this corresponds to a second order phase transi- 
tion, where the correlation length in lattice units diverges. The basic concepts and 
definition of lattice gauge theory are reported in chapter 1 . 

An important advantage of the lattice formulation is that the path integrals which 
correspond to expectation values of physical observables can be computed numeri- 
cally via Monte Carlo simulations. The idea is to generate samples consisting of a 
large number of field configurations according to the Boltzmann distribution and to 
evaluate the observables as sample average. 

Since the first works by Creutz [12] on SU(2) this method was successfully applied 
to lattice pure gauge theories and in recent years to full QCD with two and three 
flavors of dynamical fermions; in this period, increasing precision and reliability due 
to theoretical improvement of the techniques has been accompanied by an enhanced 
computing power. 

Due to the non-perturbative nature of the formulation, lattice field theory pro- 
vides a very powerful tool to be adopted for the investigation of low energy features 
of QCD; the final goal is to test the Standard Model by evaluating from first prin- 
ciples physical quantities which one may compare to the experimental results. 
Within this challenging task, lattice gauge theories can give important indications 
regarding basic questions related to strong interactions. 

For example, it is expected that at high energies perturbation theory gives reliable 
estimates of physical quantities, but how high the energy must be is still subject of 

debate. The first evident problem of the pcrturbativc expansion for QCD is that 
it is described by an asymptotic series: this means that after a certain order the 
perturbative expansion looses its predictive power. 

Apart from this, a second problem is that the accuracy of the prediction is highly 
dependent on which renormalization scheme is adopted to define the observables. It 
turns out that different ways to subtract the divergences lead to different definitions 
of the running coupling in terms of which the perturbative series can behave quite 
differently. 

The question whether perturbation theory is really applicable can be relevant in 
many practical cases; for example, the study of the perturbative evolution of deep 
inelastic scattering structure functions starting at a renormalization point below 
IGeV was strongly debated. Another popular example is the experimental mea- 
surement of the strong coupling ct^ itself [13], which needs perturbation theory as 
theoretical input. The main phenomenological test of the applicability of pertur- 
bation theory is to check whether the determinations from different processes are 
compatible. From the experiments it is indeed difficult to study a process system- 
atically as function of the energy. 

What lattice gauge theory can do in this context is to furnish a non perturbative 
evaluation of the running coupling in order to compare the results with the predic- 
tions of perturbation theory and to establish at which energy one is allowed to apply 
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Here one has to overcome a stringent limitation: the lattice spacing a can not be 
made arbitrarily small and on the other hand the physical extent of the system must 
be large (> 1.5fm) in order to approximate the condition of infinite volume. With 
this limitation, one currently reaches a^^ ~ 2GeV for Nf = 2 flavors of dynamical 
quarks. A possible alternative has been proposed in [14], where one considers a finite 
size effect as the physical observable which defines a renormalized coupling. This 
method has been applied successfully in the framework of the Schrodinger func- 
tional scheme for the evaluation of the running coupling for Nf — [15-17] and 



An important result was also the computation of the so-called Lambda-parameter 
in the theory with iVf = [15]: this quantity plays an important phenomenological 
role since it fixes the energy scale of QCD. 

For this particular scheme the scale dependence of the running coupling has been 
found to be in good agreement with perturbation theory for a < 0.3. These studies 
showed also that the perturbativc scaling regime and the low-energy domain of the 
theory are smoothly connected, with no complicated transition region. 

It is important to check these properties in many independent ways by comparing 
several quantities computed non-perturbatively with their perturbative prediction. 
An interesting quantity for this purpose is the static quark potential, which defines 
the energy of gauge fields in presence of two static color sources separated by a dis- 
tance r and has been studied very intensively in the last years. 
However, although many non-perturbative computations of this observables have 
been performed, there was up to now no convincing investigation in the short dis- 
tance regime. With this goal we decided to evaluate the static quark potential in 
the Nf = case (for our purposes this means the pure Yang- Mills theory), for which 
it is possible to simulate large lattices; this allowed the investigation of the small 
distance region and the continuum extrapolation of the results obtained at finite 
lattice spacing. 

At high energies, corresponding to small distances between the static quark-antiquark 
pair, the asymptotic freedom justifies the application of perturbative methods. Af- 
ter subtracting a constant part, the perturbative potential can be written in the form 



which defines the running coupling in the so-called V scheme. 

This definition is restricted to perturbation theory since from the non-perturbative 
point of view the subtraction of the self-energy is not well defined. 
Alternatively, one can define the coupling in the qq scheme through the force 



Nf = 2 [18,19]. 



V{r) 



3 r 



r 



1 



(2) 



F{r) 



dV _4: Q;qq(/i) 




1 



(3) 



dr 3 



In perturbation theory, the potential is known at two loops [20-22]; however, the 
reliability of perturbation theory has been doubted at distances as short as O.lfm, 
and the presence of large non-perturbative terms has been argued [23,24]. 
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We started our investigation by constructing in the continuum a perturbative pre- 
diction for ayil-i) and Q;qq(/i). By using tlie (non-perturbative) knowledge of tlie 
Lambda-parameter in units of the low energy scale ro ~ 0.5fm [25], it is possible to 
extract a parameter-free prediction for the running coupling via the renormalization 
group equation (RGE). The only ingredients that one needs are, for a given renor- 
malization scheme, the corresponding Lambda-parameter and the coefficients of the 
/3-f unctions. The Lambda-parameters in different schemes are connected exacly by 
one-loop relations and the /3-function coefficients can be taken from the literature. 
With these ingredients we evaluated the 2-and 3-loops RG perturbative predictions 
for ayi^fi) and aqq{fi). 

We want to stress that we are dealing with a single-scale problem and hence we ex- 
pect that this is the best perturbative prediction, since there are no free parameters 
related to the energy scale to be fixed. 

By analyzing the difference between the 2-and 3-loops behavior one can argue that 
for the V scheme the perturbative expansion appears to be applicable only up to 
a ~ 0.15, while for the qq scheme the series is well behaved until a ~ 0.3. This 
analysis gives already important indications on the quality of the perturbative pre- 
diction within these two different schemes. To clarify whether unexpectedly large 
non-perturbative terms are present in these observables was the strongest motivation 
of our computation. 

First of all, we tried to give a reasonable phenomenological definition of a "large non- 
perturbative term". We then evaluated the potential and the force on the lattice in 
the region 0.05fm < r < 0.8fm. 

The first step was to set the scale by specifying a dimensionful low energy observable 
in terms of which all other dimensionful quantities can be expressed. This observable 
must be chosen very carefully and must be evaluable with good precision. 
In particular the length scale ro ~ 0.5fm [25], which is defined through the force, 
turned out to be a good choice. Due to the fact that we are interested in the short 
distance region, we decided to introduce a smaller reference scale Tc ~ 0.26fm also. 
Once renormalized quantities have been constructed, we performed the continuum 
extrapolation of the potential and the force; in this procedure we were also able to 
test the theoretical expectation on the leading lattice spacing errors to F(r). 
Finally we compared our continuum results with the perturbative predictions and 
we discussed the domain of applicability of perturbation theory for the different 
renormalization schemes. 

As already mentioned, while perturbation theory describes the small distance 
regime, the long range properties are of non-perturbative nature. An intuitive picture 
of the confinement mechanism is given by an effective bosonic string theory [26] ; a 
tube of color electric fiux is formed between the quark and the antiquark and in the 
limit of large separation r one expects that the fiux tube behaves like a string with 
fixed ends. 



With this assumption, the static quark potential has the asymptotic expansion for 
large r [27, 28] 
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where 7 is an universal factor depending only on the dimensions of the system and 
a defines the string tension. 

This expression furnishes another parameter-free prediction of the potential and also 
in this case the non-pcrturbative results obtained via lattice calculations can be use- 
ful to validate this picture. 

The problem here is that it is very difficult to evaluate the static quark potential 
at large distances with a good statistical precision; special numerical techniques are 
required and this was beyond the purpose of our work. 

We compared however our results with eq. (4) also including a higher order term in 

the effective theory, which has been estimated in [29]. 

Chapter 2 is devoted to the static quark potential, the details of the lattice eval- 
uation and the comparison with the predictions of perturbation theory and string 
model. 

The second part of this work is dedicated to the study of the properties of dif- 
ferent formulations of the gauge action on the lattice. 

The simplest way to formulate a lattice gauge action is the one introduced by Wilson: 
one sums over plaquette terms, which are gauge invariant products of link variables 
along the smallest closed loop that one can build on a cubic lattice. 
Within this framework one can show that the leading discretization errors that one 

introduces are of order O(a^). 

Due to the practical limitations which prevent to study systems with very small 
lattice spacings and the impossibility in many cases to perform a continuum ex- 
trapolation of the physical observables, it would be highly desirable to work with a 
formulation in which the lattice artefacts are reduced. This problem is even more 
drastic in the fermionic action, where the leading discretization errors are linear in 
the lattice spacing. 

For this reason in the last years there were big efforts in the formulation of both 
gauge and fermionic improved actions. The basic idea is to add appropriate com- 
binations of irrelevant operators to the lattice action and to tune their coefficients 
such that the discretization errors are reduced. Using the language of statistical 
mechanics, the concept of universality ensures that the influence of irrelevant op- 
erators on the critical behavior is negligible and hence one recovers the desired 
physical continuum limit. At finite lattice spacing, the irrelevant operators governe 
the discretization errors, also called scaling violations ^, associated to renormalized 
dimensionless quantities. 

Several approaches were studied with this purpose. A popular one is based on the 
Symanzik picture [30] and has the goal to eliminate the leading lattice artefacts 
proportional to a for the fermionic action and to for the gauge action. 
This procedure has been applied to the gauge sector yielding the so-called Liischer- 
Weisz action [31] and to the fermionic action by defining the so called Sheikholeslami- 
Wohlert (clover) action [32]. 

A different strategy is based on Renormalization Group (RG) considerations; a pop- 
ular example is the Fixed Point (FP) action [33], formulated both for gauge and 



^Note that in this context scaling is used with a different meaning than before and it is related 
to the behavior of physical quantities by approaching the continuum limit. 
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fermionic field. 

Other examples are the Iwasaki [34] and DBW2 [35,36] gauge actions; the philosophy 
here is to perform blocking transformations to study the evolution of the couplings 
in a restricted parameter space, with the final goal to define a theory which lies close 

to the renormalized trajectory^ where no lattice artefacts are expected to be present. 
In particular, these actions are restricted to a two-parameter space and they were 
constructed by adding rectangular (1x2) loops with a certain coefficient to the usual 
plaquette action. 

The Iwasaki action has been used for instance by the CP-PACS collaboration [37] for 
an advanced computation of light hadron spectrum and RG improved gauge actions 
have been considered to be particularly suited for next simulations on Ginsparg- 
Wilson/domain wall fermions. 

Before one starts to use extensively these alternative actions, it is important to in- 
vestigate their properties, starting by checking the fundamental one which is the 
universality. 

Moreover, the scaling behavior has to be tested for a possibly large number of ob- 
servables in order to establish how efficient is the improvement. 
We decided in particular to study the critical temperature for deconfinement, Tc, 
for DBW2 and Iwasaki action (chapter 3); by normalizing the quantity through the 
string tension T^/ ^/o", it was observed [38] that even in the continuum limit the re- 
sult obtained with Iwasaki action deviates from what was evaluated with the Wilson 
plaquette action. 

A possible reason for this discrepancy could be a violation of universality, but we 
suspected that a more natural explanation is that the string tension is difficult to 
determine, and it is preferable to use the scale Tq to reliably set the scale in pure 
Yang-Mills theory. 

In order to check whether this is the case, we evaluated Tq at the values of the critical 
couplings that can be found in the literature and built the renormalized quantity 
Tj-Q for the Iwasaki and DBW2 action. 

In our evaluation we remarked that for these actions, and more generally for ac- 
tions containing additional loop terms, the physical positivity is violated, as already 
pointed out in [39]. This fact spoils the applicabihty of the variational method, 
which in many cases is necessary to extract the physical observables. In appendix 
B the properties of the transfer matrix for improved actions are discussed and the 
presence of unphysical states is studied by investigating the location of the poles in 
the propagator at the lowest order of perturbation theory. 

Finally we considered the glueball masses as additional quantities to test the 
properties of RG gauge actions (chapter 4). 

Glueballs represent one of the most fascinating predictions of the gauge sector of 
QCD; the mass of the lightest (0"^+) glueball is moreover particularly promising to 
study the scaling behavior since several calculations with the Wilson action [40-43] 
showed large lattice artefacts. We evaluated the and 2++ glueball masses for 
the Iwasaki and DBW2 actions for several values of lattice spacings by using an ap- 
propriate basis of operators from space- like Wilson loops up to length eight. Then 
we made use of our evaluation of Tq to form mo++ro and m2++ro. 
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There are several open questions that could be addressed and that have to be dis- 
cussed in order to find an optimal combination of gauge and fermionic actions to 
obtain reliable results in future simulations of QCD. Some of them will be discussed 
at the end of chapter 4. 



Chapter 1 

Lattice gauge fields 



Although QCD is the best candidate for a theory of the strong interactions, the 
quantitative information that can be extracted is still limited. It turns out that 
perturbation theory fails to reproduce many of the essential low-energy properties. 
The lattice formulation is one of the most elegant and powerful non-perturbative 
methods, and was originally introduced to investigate the confinement features. Its 
power comes mainly from the possibility to cvahiatc physical quantities "exactly" 
via numerical simulations, yielding a method to test whether QCD provides the 
correct framework for describing strong interactions. The QCD Lagrangian has six 
unknown input parameters: the coupling constant and the masses of the up, 
down, strange, charm and bottom quarks. ^ In the lattice regularization, once these 
free parameters have been fixed for example in terms of six measured hadron masses, 
then the (Euclidean) properties of the other particles made up of these quarks and 
gluons have to agree with experiments. Thus the lattice formulation allows the 
evaluation of physical observables from first principles. However, there are practical 
limitations, mainly due to the fact that the accessible lattice volumes and resolutions 
are restricted by the available (finite) computer performance and memory. 
In this chapter we will concentrate on the gauge sector of QCD and describe different 
lattice formulations of SU(A^) gauge theory that will be used later for the extraction 
of physical observables. 

1.1 Wilson action 

The starting point for a lattice formulation of quantum field theory is the path 
integral formahsm in Euchdean space-time [45], [46]. In this context it is possible 
to establish a useful analogy between Euclidean quantum field theory on a lattice 
and statistical mechanics, which is also the basis for numerical simulations. In the 
following we will consider pure SU(A^) gauge theories on a 4 dimensional hypercubic 
lattice 

A = aZ-^ = {x\x^/a eZ}, /x = 0,1,2,3, (1.1) 

-"^The top quark is considerably heavier than the other quarks {nit = 174.3 ± S.lGcV [44]), so 
that its bound states decay too fast to be detected. For this reason the top quark is not included 
here. 
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where a is the lattice spacing. The lattice introduces an ultraviolet cutoff and 
provides the only known consistent non-perturbative regularization of non-Abelian 
gauge theories. The Fourier transformation of functions on the lattice is periodic in 
momentum space with periodicity 27r/a, therefore all momenta can be restricted to 
the first Brillouin zone 

B={p\--<p,<-], (1.2) 
I a a) 

and the momentum cutoff is ir/a. 

An SU(A^) lattice gauge field is an assignment of a matrix U{x,fi) G SU(iV) to 
every lattice bond with endpoints x and {x + afi), where fi denotes the unit vector 
in the positive /i direction. Under a gauge transformation V{x) the link variables 
transform as 

U{x,i^) -^U'{x,ii)= V{x)U{x,n)V-^{x + afi), V 1/(a;) e SU(Ar). (1.3) 

A particular gauge-invariant object one can construct on the lattice is the trace of 
the product of link variables along a closed curve. These loops can be of arbitrary 
size and shape, and can be taken to lie in any representation of SU(A'^). The simplest 
example is the plaquette W^^^^, a (1 x 1) loop 

Wl^^{x) = Tr {U{x, n)U{x + a/i, v)U~^{x + av, //)C/~^(x, v)] . (1.4) 

The action which has been proposed by Wilson [11] for pure lattice gauge theory is 
defined in terms of these plaquette variables 

^ - /^EE ji - ^^^K^\^)] ■ (1-5) 

Let A^{x) = —igQA'l^{x)Tb, where are the generators of the group, be a Lie algebra 
valued vector field defined on the lattice and let U {x, ji) — e~"-^^'^^'> be a smooth field; 
eq. (1.5) takes the form 

^ = E ci'l^F,^{x)F^'{x) + 0{a% (1.6) 

X 

where 

F^,{x) = A^A^ix) - AlAi,{x) + [Ai,{x),A,{x)] 

and Aj^f{x) = ^{f{x + aft) — f{x)) is the lattice forward derivative. The leading 
term of eq. (1.6) for small a coincides with the Yang-Mills Euclidean action 

Sym = / d'xTrF,,F,, = ^J d'xF^^F^^, (1.7) 

if we set (3 = ^ and identify go with the bare coupling constant of the lattice theory; 

thus the Wilson action gives the desired classical continuum limit. 

After having defined the field variables and the action, the next step to quantization 



1.2 Continuum limit and improvement 



11 



is to specify the functional integral. Given an observable 0[U] which is in general a 
gauge invariant function of the link variables, its expectation value is given by 



(^) = ^J D[Up[U]e- 



■S[U] 



where 



D[U]e-^^^^ (1.9) 



and D[U] = Ylx n^^i-'^^ f^) invariant group measure or Haar measure. In 

analogy to statistical mechanics, Z is called partition function. The high dimensional 
integral in eq. (1.8) can be evaluated by means of a (stochastic) Monte-Carlo method 
as an average over an ensemble of n representative gauge configurations 

(o) = i|:oK] + Ao(-L). (1.10) 

where AO is the statistical error; the method represents an exact approach in the 
sense that the latter can in principle be made arbitrarily small by increasing the 
sample size n. 



1.2 Continuum limit and improvement 

Once one has evaluated a certain observable in lattice gauge theories, it remains the 
issue to understand how this is related to the continuum physical world. In order 
to obtain the theory in the continuum, the lattice spacing has to be sent to zero; 
at the same time, the cut-off goes to infinity and one has to construct renormalized 
physical quantities which remain finite in the continuum limit. In practice one has to 
compute the quantities of interest at different values of the lattice spacing a and to 
extrapolate the results to a = 0. An obvious limitation in this procedure is the fact 
that it is not possible to perform numerical simulations at arbitrarily small lattice 
spacings. In current calculations, one reaches a ~ 0.05fm for pure gauge simulations 
and the quenched approximation, and a ~ O.lfm for full QCD. 
Then the improvement of the lattice formulation turns out to be a very important 
topic, and in the last years there were a lot of efforts in this direction, both for the 
gauge and the fermionic action. The first requirement for an improvement program 
is a detailed theoretical understanding of the approach to the continuum limit, then 
extensive numerical studies are needed to confirm or disprove the expected behav- 
ior. The properties that a "good" improved action has to satisfy are first of all the 
reduction of the lattice artifacts, the better restauration of rotational and internal 
symmetries; in addition, for the study of topological properties, it would be desir- 
able to work with a lattice action where small distance dislocations are excluded and 
only physical instantons are present. In the following sections several improvement 
procedures will be briefiy explained. 
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1.2.1 The Symanzik program and the Liischer-Weisz action 

In 1983 Symanzik [30] suggested that the lattice theory can be approximated by a 
low-energy continuum effective theory, and the associated action can be written as 

Ses = J d^x {Coix) + aCi{x) + a^C2{x) + ...} , (1.11) 

where Cq denotes the continuum Lagrangian and £fc, with k >1, are linear combi- 
nations of local operators of dimension A + k. From the list of all possible such fields 
only those which are invariant under the symmetries of the lattice theory have to be 
included. 

In the effective continuum theory, renormalized lattice fields are represented 
through effective fields of the form 

0eff = 00 + 001 + a^02 + •■■ (1-12) 

The so called Symanzik improvement programme consists of removing lattice 
artifacts organized as a power series expansion in a, by adding irrelevant terms to 
both the action and the operators. For pure gauge theory the corrections begin at 
O(a^), since there are no dimension 5 operators. It turns out that there are only 
3 gauge- invariant dimension 6 operators and, if the Symanzik approach is correct, 
these should be sufficient to cancel the O(a^) effects in the Wilson action. The 
improved action for SU(3) then takes the form 

5 = /3^5^Q(l-^ReWc|, (1-13) 



=0 C&Ti 



where C are oriented paths belonging to sets Tj, i = 0, 1, 2, 3 of structurally equiva- 
lent curves, which are indicated in Fig. 1.1 respectively by (a), (b), (c) and (d) 
Tq— set of curves enclosing one plaquette; 

Ti= set of planar curves with perimeter 6a enclosing two plaquettes; 
T2iT^— set of non-planar curves with perimeter 6a. 

Note that the first term corresponds to the usual plaquette action and the number 
of the other classes matches the number of independent operators of dimension 6. 
The coefficients Cj have to be fixed such that on-shell quantities do not have O(a^) 
corrections; at the classical level, they can be easily determined yielding the following 
results [47], [48] 

Co = ^ = l-8ci, (1.14) 

Cl = (1-15) 
C2 = 0, (1.16) 
C3 = 0. (1.17) 

In quantum field theory the coefficients are subject to radiative corrections; the 
coefficients Cj have to be determined non-perturbatively to fully eliminate the O(a^) 



1.3 Renormalization group improved actions 



13 



(a) 



(b) 



(0 (d) 
Figure 1.1: List of paths which contribute to the action eq. (1.13). 

artefacts, and this has not been achieved yet. The 1-loop improved action goes under 
the name Liischer-Weisz action [31] and has coefficients 



co(^o) 

ci(^o) 

C2(^o) 
C3(^o) 



---0.02521^0^ 



-0.00441^0, 
0, 



(1.18) 

(1.19) 

(1.20) 
(1.21) 



with leading corrections of order 0(gQa'^). 

1.3 Renormalization group improved actions 

In recent years an alternative improvement strategy based on renormahzation group 
(RG) considerations was developed. In general a quantum field theory is defined 
over a large range of scales from low (relative to the cut-off) physical scale up to 
the cut-off which goes to infinity in the continuum limit. The goal of a renormaliza- 
tion group transformation (RGT) is to reduce the number of degrees of freedom by 
integrating out the very high momenta and taking into account their effect on the 
remaining variables ^. 



This can be achieved for example introducing a blocked lattice with lattice spac- 
ing a' = 2a and blocked fields as average of the original fields; integrating out the 
latter keeping the block averages fixed, one obtains a new action S' which describes 
the interaction between the block variables. The important observation is that the 
partition function in unchanged, so that the long-distance behaviour of the Green- 
functions is expected to remain unchanged as well. On the other hand, the lattice 
spacing is increased by a factor 2 and the dimensionless correlation length, which 



^This is the Wilson approach to Renormalization Group [49] 
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represents the inverse mass of the hghtest state in the theory, is reduced by a factor 
2 

^ 2 

In momentum space, performing this transformation is equivalent to integrate out 
all momenta between a given scale q/2 and q. 

In general, the new action S' will contain every possible interaction term, even if the 
original action S had a simple form. It is then useful to start with a general action 
of the form 

S{cl>)^Y.^M). (1-22) 

a 

where 4> are the original fields ^, 0^ are all possible operators (interactions) and 
the corresponding couplings. The transformed action where x are the blocked 

fields, is also expanded in terms of these operators 

S\x) = Y.^'Ja{x)- (1-23) 

a 

A RGT can be seen as a motion in the couphng constant space: {K^} — > {K'^}. 
Under repeated RGTs a coupling constant flow (RG flow) is generated 

{K.} - {K) - {K'\} - (1.24) 
while the dimensionless correlation length is reduced at every step 



2^4' 



(1.25) 



The physical theory is preserved under each RGT, and thus along the whole flow. 
There will be special points called fixed points (FP) defined by 

{K] {K}, (1.26) 

where the theory reproduces itself under a RGT. Due to eq. (1.25), at a fixed point 
the correlation function must be or oo. 

Now consider a point {Ka} close to the FP {K^}, with AK^ — — small. 
Applying a RGT and expanding around {K*} we have 

AK^ ^ AK'^ ^K-K ^J2^apAKp + 0{{AKf), (1.27) 

n 

with 

Tap = p.jj' \k=k*- (1-28) 

Let us denote the eigenvectors and eigenvalues of the matrix T by ha and \a respec- 
tively: 

Y,%phl = \-hl, a = 1,2... (1.29) 

H 

^In this example for simplicity we can use scalar fields defined on the lattice sites. 
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The eigenvectors h"' define the eigenoperators 



(1.30) 



which define a new basis in terms of which the action can be expanded 



(1.31) 



a 



where c" arc the corresponding coefficients and are called scaling fields. By analogy, 
/i" are called scaling operators. 

Iterating the RGT, one obtains that the couphng goes over (A")"c" after n RGT 
steps. For |A"| > 1(|A"| < 1) the coupling |c"| is increasing (decreasing) under re- 
peated RG steps. The corresponding interaction /i" is called relevant (irrelevant). 
For I A" I = 1 the behavior of the operator (called marginal) is decided by higher order 
corrections in eq. (1.27). 

These general ideas can now be applied to Yang- Mills theory; due to asymptotic 
freedom, the continuum renormalizcd theory is obtained as (70 ^ 0(/9 — 00) in which 
limit the correlation length goes to infinity, the lattice spacing a — > 0, the resolu- 
tion becomes infinitely good and the cutoff artifacts disappear from the predictions. 
Introducing an infinite dimensional coupling constant space. 



under a RGT we will have {/?, f^i, •••} {P' , K[, K2, ...}. The = 00 hyper- 
plane is the so called critical surface (where ^ = 00). Starting from a critical point, a 
RGT produces another critical point. On the critical surface there is a FP with co- 
ordinates K^,K2,.... The FP defines a basin of attraction, i.e. a set of critical points 
that converge to it under RGT. The basin of attraction determines the universality 
class and all theories belonging to them show the same long-distance behavior. The 
fiow along a relevant scaling field, whose end-point is the FP, is called renormalized 
trajectory (RT). In Yang- Mills theory, all the eigenoperators lying in the critical sur- 
face are irrelevant, and there is one marginal direction, which becomes relevant by 
higher order corrections, pointing out of this surface. The RT is the attractor of all 
flows terminating in the critical points that lie in the basin of attraction of the FP. 
Along the RT there are no scaling violations as, by construction, all irrelevant fields 
are zero. Thus, simulations done using an action along the exact RT will reproduce 
the continuum physics without discretization errors. Simulations on finite lattices of 
size L are done in a region of coupling space where <^ L, that is they lie on fiows 
that terminate at critical points that may be quite far from the FP. Corrections to 
scaling along these fiows may therefore be large. A method to reduce these correc- 
tions is to adjust the action to lie on another flow that starts closer to the FP; since 
the flows are attracted by the RT, its location may be estimated by starting with 
the simple plaquette action and performing some blocking transformation. 
A usual assumption is that the FP action is local, that is the strength of the cou- 
plings fall off exponentially with the size of the loops corresponding to the operators; 
nevertheless, it still involves an inflnite number of couplings, but in practice one is 
forced to truncate the action to a flnite number of interactions. 
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1.3.1 RG gauge improved actions in 2-parameter space 

Investigations on the Renormalization Group in (T-models and Yang-Mills theories 
showed that a satisfactory treatment of RGT's can be achieved by restricting the 
action to a two-parameter space. The most popular examples are the Iwasaki and 
DBW2 actions where, in addition to the usual plaquette term, planar rectangular 
(1x2) loops are included 

^ = /^E (^oE {i - l^^K^'i^)] + ^1 E {i - l^^K^^'i^)]) ' (1-33) 

with the normalization condition cq = 1 — 8ci. Note that eq. (1.33) has the same 
form as the tree-level Symanzik action eq. (1.13), eq. (1.14). The coefficients Cq,Ci 
are now determined by RG considerations. 

Iwasaki [34] introduced a block-spin transformation with the purpose to obtain 
a gauge action that after a few blockings comes close to the RT. The strategy is 
the following: (i) one assumes to know the lattice action on the RT; (ii) the action 
can be expanded in perturbation theory ; (iii) the correlation functions on the RT 
are calculated in perturbation theory. In principle, if one would know the (iii) then 
one can follow the opposite direction (iii) (i) and determine the action on the 
RT. In practice, one truncates eq. (1.32) to a two-parameter space and chooses as 
improved action a lattice action which is located near the RT. Iwasaki obtained for 
the coefficients in eq. (1.33) the numerical value Ci = —0.331. 

The QCD-TARO collaboration has investigated the flow of the Wilson action 
under blocking [35,36]. In particular, they used a blocking transformation with scale 
factor b — 2 and determined the couplings on blocked lattices using the Schwinger- 
Dyson equation. The so called DBW2 action is obtained after a Double Blocking of 
the Wilson action in 2-dimensional coupling space; starting from the coupling {f3 = 
6.3) on a 32^ x 64 lattice, one obtains Pi^i = Pcq ~ 7.986, /9ix2 ~ Pci ~ —0.9169; 
the corresponding coefficient in eq. (1.33) is then ci = —1.4088. 

To summarize, in a two-parameter space, the coefficient Ci in eq. (1.33) takes 
different values for various choices of alternative actions 

{—1/12 Symanzik, tree level impr. 
-0.331 Iwasaki, RG (1.34) 
-1.4088 DBW2, RG 

Notice that the strength of the rectangular loops for RG improved actions is larger 
than what is expected in order to cancel the O(a^) at tree level, i.e. this method 
suggests an over-correction; this feature will be discussed later in chapter 3. 



1.3.2 FP action 

Hasenfratz and Niedermayer [33] used the framework of the renormalization group 
to propose a classically perfect action. The basic idea is to consider the FP described 
before, to evaluate the action at the FP, and then allow (5 to move away from oo. 
This does not represent a RG flow but defines an action for every value of p. This 
defines a "classical perfect action", in the sense that its classical predictions agree 
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with those in the continuum, independently of whether the lattice is fine or coarse. 
On can say that the FP action is an on-shell tree-level Symanzik improved action to 
all orders in a. The detailed form of the FP action and the RT will depend on the 
form and parametrization of the block transformation; the theoretical properties of 
the FP action and RT are, however, independent of these details. 
The next step is the observation that in asymptotically free theories like 4-dimensional 
SU(iV), the FP lies aX j3 = oo and in this limit the path integral describing the kernel 
of the blocking transformation can be calculated in the saddle point approximation. 
This leads to an equation in classical field theory which determines the FP action. 
However, this step is too demanding to be performed in a Monte Carlo simulation and 
one needs a sufficiently fast but at the same time accurate method to approximate 
the FP action. The first works on SU(3) gauge theory proposed a parametrization 
which used plaquettes, rectangular loops and their powers. Afterwards, [50] a new 
parametrization using plaquettes made by original and smeared ("fat") links was 
investigated. 
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Chapter 2 

The static quark potential 



The potential between a static quark and anti-quark is an object of considerable 
theoretical interest that has been studied for almost 30 years. It defines the energy 
of gauge fields in the presence of two static color sources separated by a distance r. 
From the phenomenological point of view, in the limit where the mass of a quark 
q becomes large compared to typical QCD scale, bound states qq are expected to 
be described by an effective nonrelativistic Schrodinger equation; the nonrelativistic 
potential is given by the static quark potential between charges at distance r. In 
the real world, this description is expected to be approximately valid for the bb and 
maybe for the cc spectra. ^ Due to asymptotic freedom of Yang-Mills theory, one ex- 
pects at short distances a Coulomb-like behavior and the reliability of perturbation 
theory. At large distances quark confinement shows up and perturbative methods 
are no longer able to describe the behavior of physical observables; in this context, 
lattice gauge theories play an important role providing a full non perturbative ap- 
proach. Lattice results can be used to investigate non-perturbative features of the 
static quark potential as well as to test the range of validity of perturbation theory, 
which is still an important phenomenological issue which has to be investigated more 
precisely. 

In this chapter the static potential will be defined and the parameter-free pre- 
dictions from perturbation theory and the bosonic string picture will be discussed. 
Then the lattice formulation, the simulation details and the continuum extrapolation 
will be explained and the results will be discussed and compared with perturbative 
and string predictions. 

2.1 Definition 

In the framework of the Euclidean functional integral, observables are gauge invariant 
functions of the field variables; as already mentioned, particular choices are the traces 

^The fact that the interaction energy within a heavy quark bound state whose effective Hamil- 

tonian contains a kinetic term will approach the static potential defined through the Wilson loop 
in the limit of infinite quark masses is still not clarified. Hence, in principle one should distinguish 
between static and heavy quark potential. 
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Figure 2.1: A Wilson loop C^,*.. 



of parallel transporters around closed loops, which are given on the lattice by 

TrjU{C), (2.1) 

where Tr^ denotes the trace in some representation j of the gauge group and U (C) 
indicates the product of link variables along the closed curve C. The expectation 
values of these loop variables 

W{C) = {TtU{C)) (2.2) 

are called Wilson loops; their special role is revealed by the fact that every observable 
which depends continuously on the link variables can be approximated arbitrarily 
well by expressions of the form 

J2 Yl c{Ci,...,Cn)TrU{Ci)...TrU{Cn). (2.3) 

n>0 Ci,...,C„ 

As a special case consider a rectangular loop Cr,t of side lengths r and t (Fig. 2.1); 
the static quark potential is related to the large t behavior of the corresponding 
Wilson loop in the fundamental representation ^ 

V{r) = -liml\nW{Cr,t), (2.4) 

so that 

W{Cr,t)r-t-^ooCe-'''^^l (2.5) 

The meaning of eq. (2.4) can be understood in the Hamiltonian formalism (see 
App. B), where states with several static external charges in different representations 
can be characterized by means of the corresponding behavior under gauge transfor- 
mation. Let ^'[^7] be a wave function depending on the gauge links at time zero, 
and fix the temporal gauge 

C/(x,0) = l. (2.6) 

The transfer matrix T = e"''^" is a bounded hermitian operator which acts on these 
functions. The only gauge transformations left which respect this gauge arc the time- 
independent ones; in this framework the gauge-invariant wave functions describe 

^The same static quark potential could also be defined through the correlation of Polyakov 
loops, which are products of time-like links at a certain x that wind around the lattice in the 
periodic time direction. 
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the physical states in the vacuum sector, while states with two color indices which 
transform under a gauge transformation eq. (1.3) according to 

^a/^n = V^,{x)V,-^\m,s[U]. (2.7) 

describe a state with a static quark at x and a static anti-quark at y. In the following 
we consider y — x-\-r — x-\- rp,, with fi = 1. 

Due to the gauge invariance of the Hamiltonian, sectors with different distributions 
of external static charges decouple completely. The states eq. (2.7) form an Hilbert 
space Hxy, and the static potential coincides with the ground state energy. The 
Wilson loop can be then interpreted as 

W{Cr,t) = (*|e-*=|*) = |(^(")|^)|V*^" -t^oo |(*^°)|*)|V*^('^\ (2.8) 

n 

where \E'''"'' is a complete set of eigenstates of the Hamiltonian and ^' is a state of 
type eq. (2.7) e Tigy- 

In the continuum, the static quark potential can be described also through 
eq. (2.4); the Wilson loop is defined in the continuum gauge theory as 

W{Cr,t) = (Trrexp(^£ dx'^A^^, (2.9) 

where V denotes the path ordering along the loop and the average is taken with 
respect to the Euchdean Yang- Mills action eq. (1.7). 



2.2 The static quark potential in perturbation the- 
ory 

At high energies, corresponding to small distances between the static quark-antiquark 
pair, asymptotic freedom justifies the application of perturbative methods. In this 
section we will consider the continuum Euclidean formulation of Yang-Mills theory. 

In a perturbative analysis one can show [20-22] that at least up to two loop, 
after a convenient gauge fixing, all contributions in eq. (2.9) containing connections 
to the spatial components of the gauge fields Ai{f, ±t/2) vanish in the limit t — > oo. 
The definition of the potential can then be reduced to 

V;ert(r) = -^lim ^log^TVTexp J d^xJ^A^^ ^ (2.10) 

where r = |r|, T is the time-ordering operator and 

j; = ig^5,MK^) - 5{x + r)]e{e/A - xl) (2.11) 

defines a static quark in x and a static antiquark in [x + r) that couple to the gauge 
field A„. 



22 



Chapter 2 The static quark potential 



From this definition, once a regularization scheme is fixed, for example dimen- 
sional regularization in d = 4 — e dimensions, one can perform the renormalization, 
subtracting the divergences by substituting the bare coupling with the renormalized 
one, g. For example, at two loop one obtains [20-22] 



^pert(r) = I -^,e^'%M + Ke^/, (2-12) 



with q = qi = 



V'(?)pert = --^aMs(/^){l + /l(5>/^)«Ms(/^) + /2(5,/^)"|ls(/^) + -} (^-^S) 

with Cp = I for SU(3), a = the coefficients /i, /2 are known and will be discussed 
later. The constant part Vgeif corresponds to the self energy of the static sources; 
it has to be fixed by imposing an additional condition on the potential, for example 

that it vanishes at a given distance r. 

From eq. (2.13) one can see that, once the self energy has been subtracted, the 
perturbative potential can be rewritten in a Coulomb-like form (the subscript is 
now omitted) 

t/(g) = -47rC^^, ^? = q\ (2.14) 

where av(/i) is a renormalized running coupling ("effective charge"), which defines 
the so called V-scheme, and is connected to q;ms(a*) t'Y the relation ^ 

«v(g) = Q;ms(/^) (1 + /^)«Ms(/^) + /2(?, /^)aMs(/^)^ + -) ■ (2.15) 

Alternatively, one can define a coupling from the potential in the coordinate space 
(y-scheme) 

Vir) = -C/^, M=i, (2.16) 



or use the force (qq-scheme) 



f M = = Cp^laM. (2.17) 



An important observation at this point is that only the coupling Q;qq(//) defined 
through the force is also non-perturbativcly defined, because the self energy term 
of the potential drops out. Non perturbatively it is not clear how to subtract the 
constant part, and hence ay and ay are restricted to perturbation theory. This is a 
first reason to consider the force as natural observable for comparison between per- 
turbation theory and non-perturbative QCD. Furthermore, although the couplings 
ay, ay, ctqq differ only by kinematics (differentiation, the Fourier transformation), it 
makes a big difference for the applicability of perturbation theory which one is cho- 
sen to represent the potential. This problem will be discussed in details in sect. 2. 6. 



^By matching perturbatively two schemes one always encounters a relative scale ambiguity: in 
this case for examples one has to specify the relation between q and /i. 
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In general, if g{iJ,) is a renormalized coupling, its running is described by the 
/3-function 

Pig) = i^-^g, (2.18) 

which has a perturbative expansion 

m r^-g^o -f {bo + hf + b2f + ...}. (2.19) 

The first two coefficients are universal in the sense that they do not depend on 
the renormalization scheme in which g is defined. For SU(3) gauge theory and Nf 
massless quarks one obtains [51], [52] 

*»=(SF("-H''" = (SFO'''-f''')- 

The three-loop coefficient of the /3-function can be expressed as 

bl = bl\M,=o + (eiATf + e2iV2 + e^Nf) (2.21) 

with Ci listed in Table 2.1; some of them could be taken directly from the fiterature 
[20,53-57], others such had to be computed by straightforward algebra. For 
comparison we included also the Schrodinger functional scheme ^. 

From eq. (2.18) follows that the behavior of g at high energies is given by 

where 

* = log(|). (2.23) 

and in an integration constant, the so-called A-parameter associated with the 
given scheme, with the dimension of a mass. Although classical non-Abelian gauge 
field theory does not contain any mass scale, the renormalization introduces this 
mass scale into the quantized theory: this fact is called dimensional transmutation. 
The solution of eq. (2.18), 



As = iiibof)-'"^'"'^-'/''^''' exp j' 



1^1 bl 

f3{x) b^x^ b^x 



(2.24) 



^In this discussion on different running couplings, we will include also a particular scheme used 
by the lattice community to solve the problem of matching low energy hadronic scales with high 
energy perturbative regimes. In the last ten years the ALPHA collaboration applied a finite size 
technique, the so called Schrodinger functional (SF) scheme, with the goal to compute the running 
coupling at short distances in units of the low energy scales of the theory [14, 16, 17]. The strategy is 
the following: one begins by formulating the theory in a finite spatial volume with linear extension L 
and periodic boundary conditions in all spatial directions. Next one introduces some renormalized 
coupling g{L)^ which depends only on the scale L and which can hence be considered a running 
coupling. This coupling is then computed over a range of L through numerical simulations of the 
lattice theory, using a recursive procedure which allows one to go from large values of L (where one 
can make contact with the non-perturbative scales) to the perturbative domain where L is small. 
At this scale asp{L) can be related analytically to other more commonly used schemes, like MS, 
which are defined in infinite volume. The method has been applied successfully also to the case 
with two massless flavors [19]. 
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ei 


62 


63 


MS 
V 
V 

qq 

SF 


5033 
18 

ItT^ - 247r2 - 2|9 _ 704^(3) 

T TC(3) + ivr" 

3^4 314 9 1 3569 704 a/q\ 
jTT - — TT +-g 3-C(3) 

-0.275(5) X (47r)3 


325/54 
W + ^C(3) 

W + ^C(3) + f TT^ 
2791 + ^^2 + ^C(3) 

0.0361(4) X (47r)3 





81 " 

32 8 2 
27 81 

-0.00175(1) X (47r)3 



Table 2.1: Coefficients a of eq. (2.21). 

relates the coupling g = g{fj,) to the A-parameter and will be used to obtain a pa- 
rameter free perturbative prediction for g. 

Eq. (2.15) can be generalized: given two running couplings as and as' defined in 
two different renormalization schemes, they can be matched in perturbation theory 

by 

as'{si2) = as{l2) + ff^s)as{l2y + fi'^s)as{fif + . . . , (2.25) 

where s is a relative scale which is a free parameter (for example s = q/n in 
eq. (2.15)). 

The one-loop coefficient in eq. (2.25) assumes the general form 

fi'Hs) = ^{ai + a2iVf) - 87r6olog(s), (2.26) 

where oi and 02 are listed in Table 2.2, with MS as reference scheme. The other 
coefficients can be evaluated by ff^'{s) — ff^"{s) + ff'^'{s). 



S' 




02 


V 


31/3 


-10/9 


V 


31/3 + 227£; 


-10/9 - 47£;/3 


qq 


-35/3 + 227E 


2/9 - 47^/3 


SF 


-1.25563(4)/(47r) 


-0.03983(2)/(47r) 



Table 2.2: Coefficients Ui of eq. (2.26) for S = MS 

Using eq. (2.22) and the universality of 60 and 61, one finds that the A- parameters 
in two schemes S and S' are related by 
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exactly. The fact that this matching involves only one-loop coefficients is due to the 
fact that it is related to the asymptotic behavior at /x — > cxo. 
For the 3-loop coefficient &f of the /3-function, we have 



bl' ^bl- bjl^ + ^ (fi'^is) - f^isr) . (2.28) 



This relation, together with Table 2.1 and Table 2.2, can be used to obtain the 
two-loop coefficient in eq. (2.25) 



fi 



S'S/ 



An 



(2.29) 



As a particular example, one can consider the evaluation of the coefficient 63^; the 
starting point is the relation 



Q;qq(//) = Oiyill) 1 + 2 



(2.30) 



which follows directly from the definition of the force. 

Eq. (2.30) defines the connection between the qq and V schemes; in particular, one 
can read off the matching coefficients 



/i = -26o47r, h = -26l(47^)^ 



(2.31) 



which are related to the case s = 1 in eq. (2.25) ^, with S = V and S' = qq. From 
eq. (2.28) one can easily obtain 



hi - bo: 



(2.32) 



which is the simplest way to obtain from 6^, which is known from the literature. 



In the perturbative matching of different schemes the choice of the relative scale 
factor s plays an important role for the quality of the perturbative prediction [53,58]. 
The only viable criterion for fixing s appears to be to demand that the coefficients 
fi{s) are small ("fastest apparent convergence"). For instance, the choice 



■So 



A5 



yields 



/i(^o) 

/2(so) 



0, 

bo 



bl' - bl 



(2.33) 

(2.34) 
(2.35) 



^Note that the values of 62 computed via eq. (2.28) must be independent on the factor s, so 
that one is allowed to consider a particular case. 
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and /2(so) can be shown to be close to the minimum of |/2(s)|- 
This relative scale ambiguity is related to the fact that we are matching two couplings 
at finite alternatively, if one knows the A parameter in a given renormalization 
scheme, then truncating in eq. (2.19) after the term bn-i and solving cq. (2.24) 
(numerically) for g at given fi (in units of A) defines the n-loop RG solution for the 
coupling. In single scale problems, such as the static potential depending only on 
the separation r, this is expected to be the best perturbative prediction. In contrast 
to the expansion eq. (2.25), the matching is now realized at // = oo and one does 
not need to fix a relative scale parameter s. 

The ALPHA collaboration evaluated Aj^ for A^f = from lattice calculations 
employing the Schrodinger functional scheme [15], extracting it at a sufficiently high 
scale jj, where the perturbative error is negligible and converting it to the MS scheme. 
In terms of the length scale Tq ~ 0.5fm, which will be defined in section 2.4.1, the 
result is 

AmTo = 0.602(48), (2.36) 

corresponding to 

AMs~ 238MeV. (2.37) 

This result, together with the ratios A^/A^jg- and the 3-loop coefficient of the f3- 
function 6f , allows to use eq. (2.24) for a parameter-free perturbative prediction for 
as from 2- and 3-loop RG. For completeness the numerical values of A^/A^jg and 
6f for A^f = are collected in Table 2.3, 











scheme S: 


qq V 


V 


SF 




exp(7£; — 35/66) exp(31/66) 


exp(31/66 + 7b) 


0.48811(1) 


bl X (47r)3 


1.6524 2.1287 


4.3353 


0.483(9) 



Table 2.3: Ratio of A-paramctcrs and 2-loop coefficient of the /3-function for various 
schemes for A^f = 0. These results follow from [20,53-57,59,60]. 

The couplings a^q and Oy from 2- and 3-loop RG are illustrated in Fig. 2.2, 
using the central value Aj^^o = 0.602 (the 8% overall uncertainty of this number 
corresponds to a common small horizontal shift of all curves in the figure). By 
analyzing the relative difference between the 2- and 3-loop behavior of a^^, one 
observes that the perturbative expansion appears quite well behaved up to distances 
r ~ 0.25fm. At r ~ 0.2fm one would expect the 3-loop curve to have an accuracy of 
about 10%. Since the force is completely equivalent to a^^ it is given with the same 
relative accuracy; the potential may be obtained by integration of the force. On the 
other hand, we note that the 3-loop coefficients 62 for V and ^-schemes are larger in 
comparison to b^ (see Table 2.3), in particular in the V scheme. As a consequence, 
the difference between the 2-loop and the 3-loop running coupling in this scheme 
is only small at very short distances and this perturbative expansion appears to be 
applicable only up to a ~ 0.15. This is also illustrated in Fig. 2.2. 

Furthermore, if one observes the numerical values of the matching coefficients 
/2(so) eq. (2.25) for A^f = in Table 2.4, one can notice that the SF-scheme is very 
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Figure 2.2: Running coupling obtained by integration of the RG with truncation of 
the /3-functions at 2- and 3-loop. We use Aj^ = 238MeV [15] 

close to the MS-scheme, the qq scheme is not very far, but the other schemes have 
quite large values of /2(so) in relation to the MS-scheme. In particular, the large 
coefficient between the MS scheme and the V scheme means that the direct expansion 
of the coordinate space potential in terms of q;j;js (or ccsf) is badly behaved, as it 
has been pointed out in [53,61]. 



fl'Hso) 










MS 


qq 


V 


V 


qq 


1.0653 








V 


1.6095 


0.5441 






V 


4.1303 


3.0650 


2.5208 




SF 


-0.271(10) 


-1.336(10) 


-1.880(10) 


-4.401(10) 



Table 2.4: 2-loop coefficients for s = sq = As'/As and A^f = 0. 



For the Schrodinger functional coupling, aspifJ^), it was observed in [15] that it 
is in remarkable qualitative agreement with perturbation theory for a < 0.3 and 
for a < 0.2 the 3-loop expression describes q;sf(a*) within better than 2%. It is an 
interesting question, whether the coupling in this scheme is a special case or whether 
this is a more "general property of the theory". For this purpose the perturbative 
prediction for the potential and the force will be compared with our non-perturbative 
computation in section 2.6. 
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2.3 Quark confinement and the bosonic string pic- 
ture 

Perturbation theory can be applied to gauge theories in the continuum as well as on 
the lattice. Several interesting aspects, however, are inaccessible to a perturbative 
treatment and require non-pcrturbative methods. One of these is the strong coupling 
expansion, which amounts to expanding in powers of the inverse coupling and can 
be used in the lattice formulation to investigate the region where qq is large. If one 
apphes this method to the Wilson loops cq. (2.5) [62] one finds that the leading term 
is 

W(Cr,t)^2e-'''' + ..., (2.38) 

where a — — log is the so-called string tension. By analyzing higher orders in 
the expansion it was shown [63] that the strong coupling expansion of W{Cr,t) has a 
finite range of convergence and within this range the Wilson loop obeys eq. (2.38), 
which is called the area law and establishes that the potential V{r) asymptotically 
rises linearly with r 

V{r) ~™ (jr. (2.39) 
This behavior is called linear quark confinement [11]. 

The strong coupling expansion has given an important contribution to an intu- 
itive picture of the confinement mechanism, where a tube of non-Abelian chromo- 
electic flux holds quarks and antiquarks together. In the limit of very large r, one 
can expect that these flux tubes behave like strings with flxed ends, and the expec- 
tation values of large Wilson loops could be predicted by an effective bosonic string 
theory [26]. In this theory, the Wilson loops are expressed in terms of functional 
integrals over all two-dimensional surfaces bounded by the loops. For r, t large, the 
conflguration with the minimal area will dominate; the deviation of the fluctuating 
surface from the minimal one is described by a bosonic vector fleld, whose action 
defines the effective string theory. Under this assumption, the static quark potential 
V{r) has the asymptotic expansion for large r [27, 28] ^ 

V{r) + + (J^^ , (2.40) 

where // is a mass constant depending on the regularization and 

l=-^{d-2) (2.41) 

is a quantum effect that is characteristic of the relativistic bosonic string. 7 is uni- 
versal in the sense that it depends only on the dimension d of the space-time and is 
common to a large class of effective string actions. The lattice evaluation of the static 
quark potential V{r) should allow in principle a test of the validity of eq. (2.40). 
However, it turns out that the statistical and systematic errors increase rapidly at 



^We assume as leading order of the effective string theory V{r) = ar + const. Eq. (2.40) hence 
represents the next-to-leading order expression. 
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large r, and thus it is very difficult to obtain reliable values for V{r) at large dis- 
tances. Liischer and Weisz proposed a new type of simulation algorithm [64] for the 
SU(A^) gauge theory that leads to an exponential reduction of the statistical errors 
in numerical calculations of Polyakov loop correlation functions. They evaluated 
the force V'{r) and the second derivative V (r) up to r «i Ifm [29] with sufficient 
precision to test the string predictions. For d = 3 their numerical results at large r 
are in full agreement with eq. (2.40) without need of further corrections. For d = 4, 
if one allows for a higher order term 

y(r)=<7r + /.-^(d-2)(l + ^) (2.42) 

one finds agreement in the range r > 0.5fm by setting b = 0.04fm. 

As pointed out by the authors, this good agreement seems surprising, because at 
such distances one can strongly doubt the physical picture of a thin fluctuating flux 
tube. 

A further important test is related to the energy spectrum of the excited string levels, 
which at next to leading order of the effective string theory are given by [65-67] 



TT 

En = crr + iJ,+ - 
r 



~id-2) + n 



(2.43) 



In particular, to this order the energy gap AE = Ei — V{r) is equal to ^. It was 
observed that the spectrum En does not seem to agree with the one obtained in the 
effective string picture [68-72], although this argument alone can not exclude the 
validity of the effective theory, since higher order corrections could be responsible 
for the large discrepancy. Moreover, one has to take into account the technical 
difficulties in the evaluation of the spectrum. More investigations are then needed to 
understand if the string picture is really able to describe the large distance properties 
including static quark confinement and the excited states. 



2.4 Evaluation of the potential on the lattice 

Since the original work by Creutz in the SU(2) gauge theory, many computations of 
the potential from lattice gauge theories have been performed in SU(3), but the short 
distance regime has not yet been investigated convincingly. The reason is that by 
"the potential" one means the potential in infinite volume; in practice one considers 
systems with finite extension L and one knows that for the force at distances up to 
r ~ 0.5fm, the deviation from infinite volume are small on an L'^ torus with L = 
1.5fm. With such a physical value of L, an investigation of distances of r ~ 0.05fm 
with control over 0(a/r) effects requires very large lattices. Setting the number 
of flavors A^f = (for our purposes this means the pure Yang- Mills theory), large 
lattices may nowadays be simulated. We studied the static potential in the range 
0.05 fm < r < 0.8 fm, employing a sequence of lattices up to 64^. This study of 
the potential from short to intermediate distances was performed using the standard 
Wilson action and extrapolating the results to the continuum. In chapter 3 we will 
consider the potential for improved actions with the purpose to investigating the 
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lattice artifacts, and we will restrict ourselves to coarser lattices leaving out the goal 
to perform a continuum extrapolation. 

2.4.1 Setting the scale in Yang-Mills lattice theories 

The first step in the evaluation of renormalized quantities from lattice QCD is to 
set the scale. In a pure gauge theory the coupling Qq is the only parameter in the 
Lagrangian, then one has to specify one physical quantity in order to renormalize 
the theory. It turns out that it is convenient to choose a dimensionful long distance 
observable, which introduces the scale into the theory: predictions of all other di- 
mensionful quantities are then expressed in units of this scale. They are well defined 
renormalized quantities and can be extrapolated to the continuum limit from results 
at finite resolution. 

It is clear that the scale should be chosen with care since it influences the precision 
of many predictions. From the confinement picture, one expects a linearly growing 
potential at large distances, with a certain slope which defines the string tension. 
One possible choice to set the scale is then 

(7= lim F(r), (2.44) 

r— »oo 

where F{r) — This limiting procedure is not easy to perform because statistical 
and systematic errors in the force increase with the distance. Moreover, if one 
believes the string picture described in sec. 2.3, one can expect that higher order 
terms proportional to b in eq. (2.42) can be neglected only at very large distances 
and this should be taken into account in the extraction of the string tension from 
the potential. 

In most cases, the string tension is extracted from lattice data by fitting the potential 
with a three-parameter formula 

V(r) = ai + aar + — , (2.45) 
r 

but to identify 02 with the string tension can be a source of large and uncontrolled 

systematic errors. 

For this reason it is desirable to use the force at intermediate distances to determine 
the scale. One can observe that both bb and cc spectra can be accurately described 
by one effective potential in the range of r 0.2 fm to r Ri 1 fm, and this suggests 
that the range where one has the best information on the force F{r) between static 
quarks is at distances of around 0.5 fm. One then calculates r(c) satisfying the 
equation [25] 

r^F{r) I _ = c , ro = r(1.65) ^ 0.5 fm. (2.46) 

When one is interested in short distance properties of the theory, it is convenient 
to choose a smaller reference length scale, because large distances on fine lattices 
are very expensive from the computational point of view. We decided to separately 
consider two regions, 0.05 fm < r < 0.3 fm and 0.2 fm < r < 0.8 fm, where the first 
requires very large lattices and small lattice spacings and the second one had been 



2.4 Evaluation of the potential on the lattice 



31 



simulated before on coarser lattices [73]. The region of overlap serves for calibration. 
We therefore introduced an additional scale 

rc = r(0.65) (2.47) 

and evaluated its relation to Tq. 



2.4.2 Definition of the force 

In the continuum, the definition eq. (2.46) is unique, but at finite lattice spacing one 
has to specify which discretization of the force is to be used. Given the potential 
V{r) in a certain direction, in general one has 

F(r') = [V{r) - V{r - a)] /a. (2.48) 

The naive choice for r' would be 

r' = r„ = r-|. (2.49) 

Instead of r„, we decided to use a distance ri chosen such that the force on the lattice 
has no deviations from the force in the continuum when evaluated at tree level. On 
the lattice, the force at tree level is given by 

.x__4c/o r cos{rki/a) -cos{{r -a)ki/a) 

3aJ_A2nr 4 sin^(^,/2) ^ ^ ^ 

and one would introduce ri such that 

The lattice artifacts remain only in the Olg^) terms; F{ri) is called a tree-level 
improved observable. 

In the SU(2) theory it has been observed that by adopting this definition the 
remaining lattice artifacts in the force are surprisingly small [25] and here we shall 
again find no evidence for them as long as r > 2a. 

Explicitly we have from eq. (2.50) 

{ATrr^y^ = - [G{r, 0, 0) - G{r - a, 0, 0)] /a, (2.52) 

where 

^^'^ - a I. W4E,i,sin^(fc./2) ^ ^ ^ 

is the (scalar) lattice propagator in 3 dimensions, for general f= {xi, X2, x^). In [74] 
a method was proposed to construct it directly in coordinate space using a recursive 
relation. We extended the method to 3-d finding the necessary invariant of the 
recursion relation eq. (C.8). Formulae are fisted in App. C. 
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2.4.3 Computation of the potential and force 

The evaluation of the potential from eq. (2.8) is technically difficult because of the 
presence of higher order excitations, so that the ground state becomes dominant 
only at large values of t, where the statistical uncertainties can be very large. In 
order to enhance the overlap with the physical ground state one can use smearing 
techniques, which consists of building operators which describe in our case static 
quarks and antiquarks in an "optimized" way. 

In particular, we adopted a smearing operator S which acts on the spatial compo- 
nents of the gauge fields via [75] 

SU{x, k) = Psu(3)] U{x, k) + aY, [U{x,j)U{x + aj, k)U\x + ak,j) (2.54) 

+U''{x - ajJ)U{x - aj,k)U{x + ak - ajj)] 

where Vs\j{3) denotes the projection into the group SU(3) and a is an adjustable 
parameter. ^. 

For different smearing levels / = 0, M — 1 one then constructs smeared spatial 
links according to 



In particular we chose 



In [73] ni was set to 



Ui{x,k) = S'''U{x,k). (2.55) 



a^-, M = 3. (2.56) 
2 



2 



m^- (^) . (2.57) 



According to this choice, n2 corresponds to what was estimated to be the optimal 
smearing in [76]. Since we want to compute the potential for shorter distances, where 
one can reach high precision, we adopted in general smaller values than these. 

Another difficulty is the fact that the Wilson loops eq. (2.8) decay exponentially 
for large t. In order to measure these values in a Monte Carlo simulation with 
statistical significance, also the variance of the matrix should decrease exponentially 
with t. To achieve this, one can apply a so-called multi-hit procedure [77] to the time- 
like links, which are substituted by their expectation values in the fixed configuration 
of the other field variables 

, rdC/(a;,0)C/(x,0)e-^[^(^'°)l , , 

where S[U {x, 0)] denotes the part of the action depending on U {x, 0). For a general 
observable, the substitution of U {x, 0) with U {x, 0) is allowed under the following 

'^For any 3x3 matrix V, ■Psu(3)^ denotes the matrix W € SU(3) which maximizes Re Tr(W^Ft) 
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restrictions: (i) the observable must depend linearly on the considered links and (ii) 
no pair of substituted links may be coupled through the action. 
For SU(2), an exact expression for eq. (2.58) can be obtained in terms of modified 
Bessel functions. For SU(3), there is no analytic expression and eq. (2.58) is evalu- 
ated numerically by means of Monte Carlo integration. ^ At fixed r, one then forms 
a M X M correlation matrix of Wilson loops 

Cimit) = ^Tr [Vi{Q, ri)V{ri, A + tQ)VX{tQ, A + tQ)v\Q, tO) } ^ (2.59) 

= Cml{t), 

where 

Vi(x,x + A) = Ui(x,l)Ui{x + ai,l)...Ui(x + (r -a)l,l), (2.60) 
V{x,x + m) = U{x,0)U{x + aO,0)...U{x+{t-a)0,0). (2.61) 

The correlation matrices were computed on field configurations generated by a hybrid 
algorithm with a mixture of A^or over- relaxation steps per heat bath step [78,79]. 
The ratio Not increases with the coupling /3, and in practice we adopted 

Nor ~ 1.5(ro/a). (2.62) 

In a Monte Carlo simulation the configurations are correlated, hence the measure- 
ments are not statistically independent. We decided to perform a measurement of 
Wilson loops every N a/ meas sweeps, where Nu/meas is chosen such that the time 
needed for the updating is about half of the time used for the measurement. 
By collecting the data in bins and applying the so-called jackknife procedure, it is 
possible to check that with this criterion the measurements are practically uncorre- 
lated. 

The simulation parameters for 5.7 < f3 < 6.4 are reported in [73], while for the new 
simulations in the range 6.57 < (3 < 6.92 they are listed in Table 2.5. We used 
L X L X L X T lattices, with T = L. The parameters (3 and L/a have been chosen 
such that 

L/ro ~ 3.3. (2.63) 

This condition is necessary because the physical results are obtained in the limit 
L — > oo, and by using a finite volume one has to be sure that finite-size effects are 
negligible. In [73] it was shown that eq. (2.63) is safe and no finite volume effect is 
observed. 

We evaluated the Wilson loops in eq. (2.59) for each value of the loop size r,t, 
with a < r < rmax and a < t < tmax- Due to the high precision that one can 
reach at short distances, only 40 to 75 measurements were sufficient to extract the 
observables, and with relatively modest computational effort we could extend the 
calculations to a ~ 0.025 fm where a 64^ lattice had to be simulated. 



^We evaluated the average using 15 Cabibbo-Marinari iterations. 



34 Chapter 2 The static quark potential 



Table 2.5: Simulation parameters. 



L/a P no rii n2 iVpr r^^/a tmax/a 

40 6.57 77 153 18 7 11 

48 6.69 53 106 22 10 12 

56 6.81 72 144 25 12 14 

64 6.92 94 188 29 12 16 



2.4.4 Analysis details 

The correlation matrices eq. (2.59) were analyzed using a variational method, which 
was proposed in [80,81]. For any r and to fixed (to = in practice), one solves a 
generalized eigenvalue problem 

{t,to), Xa>Xa+l, (2.64) 

which corresponds to finding the eigenvalues of the symmetric matrix 
C — C{to)~'^^'^C{t)C{to)~^/'^ . For the Wilson action, the positivity of the transfer 
matrix ensures that C{t) is positive definite for all t (see appendix B).^ One can 
prove that the energies Ea are related to the generalized eigenvalues by the relation 

aE^ = ln(A„(t - a, to)/Xa{t, h)) + 0(e-*^^«), (2.65) 

where AEa = mmp^alEa — Efi\. The static quark potential corresponds to the 
ground state energy, V{r) = Eq. 

By using this method one expects that the coefficients of the higher order terms in 

eq. (2.65) are suppressed and one can read off the energy values already at moderately 

large t. Solving the generalized eigenvalue problem eq. (2.64) is indeed equivalent to 

looking for a linear combination of operators appearing in the correlation function 

such that the state with the lowest energy has the maximal weight. 

We decided to extract the potential V{r) projecting the Wilson loop correlation 

function 

W{t) = vlC{t)vo, (2.66) 

where Vq is the ground state eigenvector computed for t = to + a in eq. (2.64) 
corresponding to the ground state, and then constructing the effective potential 

Beside the potential, one is also interested in the energy gap between the ground 
state and the first excited state 

AEo ^Ei- Eq, (2.68) 

^For gauge actions of kind eq. (1.33) containing not only the plaquette, but also other operators, 
this property is not satisfied, and the implications for the applicability of the variational method 
are discussed in chapter. 3. 



(2.67) 
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Figure 2.3: The potential evaluated at r ~ ro/2 for different couplings as a function 
of exp(— fmin^); train will be determined such that exp(— fmj„A) < 0.3. Similar 
results hold for other values of r. 



which dominates the finite time corrections in the extraction of V{r). We estimated 
them by 

a^Eo = ln(Ai(i - a, io)/Ai(i, to)) - H^t - a, io)/Ao(i, ^o))- (2.69) 

To demonstrate that these corrections are small, we computed aV{r-) from eq. (2.67) 
for r fa ro/2 and r f=:i ro/4 and we plotted it as a function of exp(— tAE'o); as 
in [73] we found that for t > tj^in excited states contaminations are very small when 
exp{—tmin^Eo) < 0.3 (Fig. 2.3), and we used this criterion to establish at which t 
the systematic errors in eq. (2.67) are negligible with respect to the statistical un- 
certainties. In [73] a continuum extrapolation of the quantity TqAEq was performed, 
yielding 

roAEo\r=ro = 3.3(1). (2.70) 

Since the gap grows towards small r from 3/ro at r = ro to f=:i 7/ro at r = O.lro, 
the short distance region is easiest in this respect. 

All errors were computed by jackknife binning, and the results for the potential and 
force are listed in appendix A, in particular in tables A. 4, A. 2. 

2.4.5 Interpolation 

Solving eq. (2.46) and evaluating the force and potential at distances r = xr^ for 
given X, requires the interpolation of those quantities. While this can in principle 
be done in many ways, it is advantageous to use physics motivated interpolation 
formulae. 
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Figure 2.4: The ratio Tc/ro for 5.95 < (3 < 6.57 (circles) including the continuum 
extrapolation (solid point). The data for 5.95 < /3 < 6.57 are taken from [73]. 



For the force we followed [25] choosing the interpolation function 

F{r) = /i + /2r-^ (2.71) 

between two neighboring points. The systematic error arising from the interpolation 
was estimated by adding a term /ar^^ and taking a third point. The difference 
between the two interpolations was added (linearly) to the statistical uncertainty. 
We observed that at least for r > OAtq the interpolation error is smaller than the 
statistical uncertainty. For small distances the systematic error increases and can 
be of the order of the statistical one or even bigger. 

The potential is interpolated by the corresponding ansatz 

Vi{r) = vi + V2r + v^r''^ . (2.72) 

Here two points are chosen such that the desired value of r is in between. For the 
choice of a neighboring third point one has two possibilities, leading to two results. 
Their difference was taken as the interpolation error. Also in this case the systematic 
errors are larger than the statistical ones at short distances while at large distances 

the situation is reversed. 

It is important to stress that the interpolation formulas eq. (2.71), eq. (2.72) provide 
excellent approximations of the r dependence only locally, and to use instead a global 
fit in a large range could be a source of large systematic errors. 

2.4.6 The ratio rc/rg and parametrization of ro/a 

Following the procedure described in the previous sections, we evaluated Vc/a and 
ro/a in the region 5.95 < /3 < 6.57 where both quantities are accessible. Figure 2.4 
shows the ratio Tc/ro for several lattice spacings as a function of (a/ro)^. No de- 
pendence on the resolution is seen within the errors of below 1%. A continuum 
extrapolation gives 



rjro = 0.5133(24) 



(2.73) 
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Table 2.6: Results for r^ja [73] and r^ja. For comparison also the values from the 
interpolation formula eq. (2.79) are shown. 



/3 


ro/a 


eq. (2.79) 


/3 


Tc/a 


eq. (2.79) ,eq. (2.73) 


5.7 


2.922(9) 


2.938 


6.57 


6.25(4) 


6.22 


5.8 


3.673(5) 


3.665 


6.69 


7.29(5) 


7.25 


5.95 


4.898(12) 


4.912 


6.81 


8.49(5) 


8.48 


6.07 


6.033(17) 


6.038 


6.92 


9.82(6) 


9.86 


6.2 


7.380(26) 


7.383 








6.4 


9.74(5) 


9.741 









and wc note that it is safe to use eq. (2.73) also at finite lattice spacings starting 
around /3 = 6.4. 

The direct determination of ro/a for 5.7 < f3 < 6.4 [73] and the new computations 
of Tc/a in the range 6.57 < P < 6.92 may be combined with rc/ro=0. 5133(24) to 
obtain ro/a in the whole range 5.7 < P < 6.92. The results are reported in Table 2.6. 
Once ro/a is known for several P values, it is useful to attempt a phenomenological 
representation in terms of an interpolating formula giving ro/a in the whole range 
of couplings. One can start from the renormalization group equation for the bare 
coupling; on the lattice the P function associated to go is defined as, 

P(9o) = -a^ = -bogi - hgl + ... (2.74) 
Solving eq. (2.74) yields the A-parameter associated to the lattice 

Klat = lim -e-V(2''o^^g)(^^^2^-fi/(2(.i)_ (2.75) 

In pure gauge theory every physical quantity with dimensions of a mass must be 
proportional to Klat in the continuum limit; therefore the mass (or a/ro in our case) 
in lattice units, which is the quantity that one evaluates in numerical simulations, 
varies as a function of go according to 

^ = ^e-^/(2''°^§)(6oyo')"''^^''°^e-'=i'°"''^°"-, (2.76) 

^0 

where A is a. constant and the terms with ci, C2, ... are due to higher order contribu- 
tions to the P function. This behavior is also called asymptotic scaling. 
The leading behavior is then given by 

^ oc e-/^/(^^^o), P^Q/gl (2.77) 

^°In [73] the simulation for /3 = 6.57 was performed with a single smearing level, yielding less 

control over excited state contaminations. As a result the error in rg/a seems to be somewhat 
underestimated. Although the statistical significance of this small effect is not clear, we decided to 
use our new data for /? = 6.57 instead. 
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Figure 2.5: Interpolation of ro/a. 



then one is justified to attempt a phenomenological representation of (ln(a/ro)) as a 
polynomial in f3. In order to avoid cancellations in the fit parameters, it is convenient 
to shift the value of j3 to {f3 — x), such that {f3 — x) is small in the interval of couplings 
that we considered. We then chose the following ansatz: 

p 

ln(a/ro) = ^afc(/3-6)^ (2.78) 

k=0 

and found that 

ln(a/ro) = -1.6804 - 1.7331(/9 - 6) + 0.7849(/3 - 6)^ - 0.4428(/3 - 6)^ 

for 5.7 < (3 < 6.92, (2.79) 

is an excellent approximation to the MC results, as shown in Fig. 2.5. The accuracy 
of ro/a in eq. (2.79) is about 0.5% at low P decreasing to 1% a,t P — 6.92. 

2.5 Continuum force and potential 

One can now determine the continuum force and potential by extrapolation of the 
MC-results at finite values of the lattice spacing to the continuum, which corre- 
sponds to a — {P — oo). 

As mentioned in sect. 2.4.5, the lattice results are first interpolated in r, then renor- 
malized dimensionless quantities are evaluated. If the lattice spacing is small enough, 
one expects renormalized dimensionless quantities to be independent of /? (and thus 
of a) ; this feature is called scaling. Any deviation from this scaling behavior is due 
to the presence of lattice artifacts. The size of the scaling violating corrections will 
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in general depend on which quantities are being considered. 

In order to approach the continuum hmit, one needs a theoretical understanding 
of the lattice artifacts; the standard framework for this discussion is Symanzik's 
effective theory, which is expected to give the asymptotic expansion of suitable lat- 
tice observables in integer powers of the lattice spacings (up to at least logarithmic 
further suppressions as a — 0) in an asymptotically free theory [30,82]. In the 
2-dimensional 0(3) a model it predicts the expansion to start at order but unex- 
pectedly numerical results are described better by a dominant linear term in a for a 
range of lattice spacing [83] . While there is no evident contradiction with the result 
of an analysis a la Symanzik which is supposed to describe the asymptofAc behavior, 
the standard picture should be tested in 4-d gauge theories and QCD as much as 
possible. 

The leading artefacts in potential differences are expected to be O(a^), but this is 
not obvious since the potential is usually defined in terms of a Wilson loop, which 
does not fall into the category of correlation functions of local fields discussed by 
Symanzik. 

In [84] a clean argument supporting this expectation is given: in short, the mo- 
tivation is that the heavy quark effective theory [85] formulated on the lattice is 
0(a)-improved without adding any additional operators to the Lagrangian [86]. The 
potential is an energy of the effective theory and thus has no linear a-effects. Of 
course, for this statement to be meaningful, the self energy has to be eliminated by 
considering potential differences or the force. 

Figure 2.6 shows the dependence of the force in units of Tc on the resolution for 
some selected values of the separation r. Fitting 

rlF{r) = rlF{r) [l + s{r/r^) x (a/re)'] (2.80) 

for fixed r/vc and for r/a > 2, the slope s is statistically not significant throughout 
our range of r and (3. Our statistical precision allows to quote 

\s{r/rc)\ < 1 for 0.4 < r/r^ < 1 . (2.81) 

Of course this 1-a bound is valid only when all details are as discussed above, in 
particular eqs. (2.48,2.52) are used to define the force at finite a. If the naive form 
r' = r — I is employed instead, the corresponding slopes s become rather large, as 
can be seen in the figure. 

A similar statement (r/a > 2), 

r2F(r) = r2F(r) I^^Jl + s(r/ro) X (a/ro)2] , (2.82) 
\s{r/ro)\ < 1.2 for 0.5 < r/ro < 1.5, (2.83) 

can be made for larger r. 

One can conclude that using a tree-level improved definition of the force, lattice 
spacing effects are below the level of our small statistical errors when one restricts 
oneself to r > 2a and to the lattice spacings considered here. By contrast, with a 

^^In general, the discretization errors parametrized by s depend also on the orientation r/r, 
which in our calculation was fixed to r/r = (1, 0, 0), however. 
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Figure 2.6: Continuum extrapolation of r^F{xrc), for x = 0.4, 0.5, 0.9 from top 
to bottom and of r'^F{xrQ), for x = 0.5, 0.6, 1.5 from top to bottom. The data 
are from our new computations and from [73]. Filled circles correspond to the naive 
value r' = r — I instead of eq. (2.52). 



naive midpoint rule for the force, the a-effects are quite sizeable. The difference is by 
construction a pure power series in a, starting with a^. For r > 2a this series is well 
approximated by the leading term. Thus both data sets in Fig. 2.6 contain the same 
essential information: full compatibility with Symanzik's theory of discretization 
errors. 

The continuum force is plotted in Fig. 2.7 using eq. (2.73) to combine the two 
regimes of r. Some data at finite (3 (in particular /3 = 6.4, 6.92) are included in the 
figure. In these cases we used our "bounds" on s to estimate that the discretization 
errors are smaller than the statistical ones. One can observe that for large distances 
the force reaches a constant value, as sign of the confinement; a detailed discussion 
on the comparison with the predictions of the perturbative theory and the bosonic 
string model will be presented in section 2.7. 

The potential contains the same physical information as the force but statistical 
and systematic errors in the lattice determination are different. From the non- 
pcrturbative point of view, the potential is defined only up to a constant, which 
represents the self energy of the static quark-antiquark pair. To eliminate this self 
energy contribution we consider Vi{r) — Vi{rc) , where we apply the same improvement 
as for the force, 

Vi{di) = V{d) , (ATTdi)-^ = G{d,0,0). (2.84) 

In analogy to ri, di is defined such that the potential has no lattice artifacts when 
evaluated at three level. The continuum extrapolation is performed exactly as for 
the force. In this case the slopes, s(r/rc), defined as above are statistically significant 
for both our smallest distances (r 0.3rc) and the largest one (r 1.5ro). As a 
consequence, at small r the combination [Vi{r) — Vi{rc)]rc has about 1% discretiza- 



2.5 Continuum force and potential 



41 



15 



10 



. continuum limit 

. 13 = 6.92 
□ i8 = 6.4 



0.5 1 
r/ro 



1.5 



Figure 2.7: The force in the continuum hmit and for finite resolution, where the 
discretization errors are estimated to be smaller than the statistical errors. 




Figure 2.8: Continuum extrapolation of Vi{r), for r/ro = 1.5, 0.6, 0.5 on the left 
hand side and for r/vc = 0.9, 0.4, 0.3 on the right hand side. 
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tion errors at j3 = 6.92 while in the large distance region these errors go up to 1% in 
[Vi{r) — Vi(rc)]ro at /3 = 6.4. These are the /9-values corresponding to the smallest 
lattice spacings available in the two different regions. 

One concludes that here, small a-cffects could also be observed for the tree level 
improved definition. As for the force -but not shown in the figure- the standard def- 
inition without tree-level improvement shows quite large a-effects compatible with 
an a^-behavior, and this confirms again the validity of the Symanzik picture. 

The continuum potential is plotted in Fig. 2.9, with some data at finite j3. In 
this figure the statistical as well as the discretization errors are below the size of the 
symbols. The results for continuum force and potential are listed in Table 2.7. 

Table 2.7: Potential and force after continuum extrapolation. 





ro-F(r) 




r/ro 


rgF(r) 


ro(Vi(r) - Vl(re)) 


0.3 




-1.676(16) 


0.5 


2.518(16) 


-0.0336(61) 


0.4 


7.30(11) 


-1.2125(93) 


0.6 


2.173(15) 


0.1989(76) 


0.5 


5.363(78) 


-0.8926(74) 


0.7 


1.951(14) 


0.4051(89) 


0.6 


4.244(53) 


-0.6475(54) 


0.8 


1.812(11) 


0.5930(99) 


0.7 


3.538(48) 


-0.4494(35) 


0.9 


1.722(11) 


0.769(11) 


0.8 


3.060(38) 


-0.27950(15) 


1.1 


1.592(10) 


1.101(11) 


0.9 


2.713(30) 


-0.13259(67) 


1.2 


1.559(13) 


1.258(12) 








1.3 


1.537(18) 


1.413(13) 








1.4 


1.537(24) 


1.567(14) 








1.5 


1.507(27) 


1.722(24) 



2.6 Comparison with perturbation theory 

The perturbative expansion in the coupling as is the most important theoretical 
tool for analyzing strong interactions effects in high energy scattering experiments. 
In practical applications however, in particular in the determination of the running 
coupling itself, it is tempting to apply the perturbative series already where the 
coupling is not so small to justify the applicability of perturbation theory. Popular 
examples are the determinations for as from r-lepton decays and the perturbative 
evolution of deep inelastic scattering structure functions starting at a renormalization 
point below 1 GeV. In order to establish the reliability of perturbation theory, it 
would be very desirable to study such processes systematically as functions of the 
energy, but either they involve a fixed energy (r-decays, hadronic Z-decays) or the 
precision is not sufficient over a larger energy range (deep inelastic scattering, e'^e~ 
total cross section, Adler function). The main phenomenological test of perturbation 
theory, therefore, is the overall consistency of the determination of as from different 
processes (see [13]). 

Complementary information may be obtained from observables that can be com- 
puted via lattice QCD, like the static quark potential. As already mentioned in 
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Figure 2.9: The static potential in the continuum limit. 



section 2.2, in perturbation theory the potential has been computed at two loops, 
but at the same time the usefulness of perturbation theory even at distances as 
short as 0.1 fm has been questioned [53,87]. Comparison between the perturbative 
prediction and the non-perturbative results will show that this is a question of a 
suitable rcnormalization scheme. When the most natural scheme, the qq defined in 
terms of the force eq. (2.17) is adopted, perturbation theory is well behaved at such 
distances. This is a rather stringent test of perturbation theory, since we can use 
the knowledge of the A-parameter together with the RG perturbative evaluation of 
the coupling, which predicts not only the scale dependence of the coupling, but also 
its absolute value. 

In an exploratory investigation, concentrated on the question, whether there are 
"large non-perturbative" terms in the potential at short distances [23] they had been 
argued to exist [24]. 

In particular, non-perturbative ultraviolet contributions to the running coupling 
proportional to A'^/q'^ were investigated. In the case of the static quark potential, 
such a contribution to the running coupling results in a term proportional to the 
quark separation, r. 

The first important step to start our analysis is to give a definition of "large non- 
perturbative term". We assume that the following is understood by this statement. 

1. A certain quantity, here the potential V{r), is considered in a region where 
its perturbative expansion looks well behaved, i.e. the n-loop contribution is 
a small correction and significantly smaller than the (n — l)-loop contribution 
(unless the latter is accidentally small itself). 

2. The difference between the full non-perturbative observable and the truncated 
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perturbative series is much larger that the last term in the series. 

With such a definition, necessarily somewhat phenomenological in character, we 
shall demonstrate that there are definitely no large non-perturbative terms in the 
potential. To the contrary, perturbation theory works remarkably well where crite- 
rion 1. is satisfied. 

In section 2.2 we established that for a^q this condition is valid up to distances 
r ~ 0.25 fm, corresponding to «qq < 0.3. In Fig. 2.10 we compare the perturbation 
theory predictions for aqq(r) with our non perturbative results, which were obtained 
from the force extrapolated to the continuum limit by the simple relation 

«qqW = lF{ry. (2.85) 

We included also additional points at finite /3, where the discretization errors were 
estimated to be smaller than the statistical ones. The 3-loop RG expression with 
^MS upper end of the error bar of eq. (2.36) is in very close agreement with 

the non-perturbative coupling. In fact the agreement extends up to values of ctqq 
where perturbation theory is not to be trusted a priori. For a^q < 0.3 criterion 1. is 
satisfied and there is no evidence for non-perturbative terms in this region. 




0.2 0.4 0.6 

Figure 2.10: Running coupling in the qq scheme. Results for the continuum limit as 
well as additional points at finite /3, corresponding to finite lattice spacing are shown. 
The perturbative curves use A^jg-ro from [15] with the dotted lines corresponding to 

the 1-0" uncertainties of this number. 

One can perform an additional test by considering the relation between as^it^) 
and Q!qq(/i) at finite /j,: using eq. (2.25) and the coefficient in Table 2.4, we define 
the 3-loop expression 

/i(asF(/i)) = asF(/u) + 1.336 [asF(/i)]^ (2.86) 

^^Note that the term proportional to q;sf(m)^ in not present due to our specific choice of the 
relative scale factor s, eq. (2.33). 
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such that 

ttqq(/i) = /i(asFWso)) + 0{[asFifJ'/so)f). (2.87) 

Then the non-perturbative values of «qq(/i) and /i(q;sf(/^/ Sq)) are compared in Fig. 2.11 
If the higher order terms in eq. (2.87) are neghgible, the two different quantities 
should agree. At a^qdj.) fa 0.3 a difference is visible but this is only about 3 x a^, 
not so far from the expected size of the next order term in the series. At Q;qq(/i) ~ 0.2, 
the difference aqq(/Li) — h{asp{fi/ sq)) is not significant at all. Wc conclude that also 
in eq. (2.87) a large non-perturbative term at short distances is excluded. 

One has to note here that the static quark potential defined through the Wil- 
son loop suffers from infrared (IR) divergences when computed at finite orders of 
perturbation theory [21]; this is related to the non-Abelian nature of SU(3). which 
allows massless gluons to self-interact at arbitrarily small energy scales. The leading 
IR singularities in the Wilson loop were argued to exist at order in [22] . These 
singularities can be regulated by a resummation, and the result is that the next order 
correction in eq. (2.87) is formally enhanced by a logarithm of a and reads [21,22,87] 

{A\og{a) + B)a\ (2.88) 

While A = ^ has been calculated [88,89], B is not known. The AQ;''log(a;) term 
by itself constitutes a small correction in the Fig. 2.11, which would slightly enlarge 
the difference between /i(Q;sF(/^/so)) and Q;qq(/x). 

Finally, Fig. 2.12 shows the static potential itself compared to different perturba- 
tive approximations. In this case we can not use the coupling ayi^fJ^) for the compar- 
ison, because it does not have a non-perturbative definition, as already mentioned 
in Sect. 2.2. The potential is defined up to a constant, which has to be eliminated 
for example by considering the quantity 

[V{r) - V{r,)]ro, 

where Tc is related to tq by eq. (2.73). 

Full line and short dashes are given by integrating the force 

V{r) = V{Q.3r,) + / dyF{y) , F(r) = Cpr^aq^il/r) (2.89) 

with the 3-loop and 2-loop RG-solution for aqg. The free constant was fixed such 
that at r = 0.3rc the perturbative prediction coincides with the non-perturbative 
result. This comparison is meaningful only with the assumption that the pertur- 
bative expansion is well behaved at that distance. Since we know that the 3-loop 
RG-solution for is accurate, this also hold for V{r) computed through eq. (2.89). 
Again, the full line moves very close to the data points (r < 0.5ro), when A^jgro at 
the upper end of the error bar is inserted. 

Another possibility is to compare the potential directly with eq. (2.16) by using 
the RG-solution for cty. From the considerations of section 2.2 we expect that per- 
turbation theory is applicable only up to ay 0.15; by solving the RG equation in 
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the V scheme with the central value eq. (2.36), one observes that at r = 0.3rc one 
has av(21oop) ~ 0.28, av(31oop) ~ 0.9. This means that one is already at those 
small distances outside the region where perturbation theory can be trusted. The 
long-dashed line in Fig. 2.12 represents eq. (2.16) with the 2-loop RG solution for 
ay; we decided not to use the 3-loop prediction due to the fact that the coupling is 
very large at the distance r ~ 0.15ro. As it was to be expected due to the missing 
stability of this perturbative expression, it fails in describing the non-perturbative 
potential. 

A similarly bad perturbative expression is the direct expansion of the potential in 
terms of Oms! ^'^^ ^^^^ comparison we used the matching equation eq. (2.25), with 
S' = V , S = MS. For the same reason explained before, we decided to use the 
2-loop RG prediction for ctj^. The dotted line corresponds to the choice s = 1 in 
eq. (2.25), while the dashed-dotted line corresponds to s — Sq (eq. (2.33)). 



I I I ' ' ' I ' ' ' I 
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Figure 2.11: Test of eq. (2.25). The uncertainty in the combination iitq has been 
translated into an uncertainty for /i(asF(/^/so)) and Q!qq(/i). The non-perturbative 
values for qsf(m) constructed from the data of [15]. Errors are smaller than the 
size of the symbols. 

In summary, care has to be taken which perturbative scheme is adopted to de- 
scribe the potential. However, perturbation theory does its best in the following 
sense. As usual in an asymptotic expansion, one should first investigate the appar- 
ent "convergence" by comparing subsequent orders and checking that they decrease 
significantly. If this is not the case, one is obviously outside the domain of applica- 
bility of perturbation theory or has chosen a bad truncation (scheme). According 
to this criterion the /9-function in the qq-scheme may be trusted up to ctqq ^ 0.3. 
Other truncations of perturbation theory for the potential are applicable for much 
smaller values of the coupling only; for example, a scheme with a large 3-loop coef- 
ficient such as V is of no use in the region a > 0.15. Therefore, perturbation theory 
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Figure 2.12: The potential compared to different perturbative expressions. Vc is 
related to ro by eq. (2.73). 



suggests that the qq scheme should be used in order to obtain a reliable perturbative 
expression. The comparison with non-perturbative results, obtained in the contin- 
uum limit of lattice simulations (iVf = 0), does confirm that such a perturbative 
analysis is a good guidehne -at least in the case at hand. At rather high values 
of the coupling, a ~ 0.3, we only expect that the 3-loop perturbative prediction 
is good to within about 10% and indeed in Fig. 2.11 one sees explicitly that that 
the truncated perturbative series has errors of this order. Similar results have been 
found for the coupling obtained with A'f = 2 in the SF-scheme [19,90]. 
If one considers the full QCD, the first observation is that, compared to Nf — 0, 
the relevant perturbative coefficients, 62 and /2(so) which are listed in Table 2.1 and 
Table 2.2, are roughly a factor two smaller in magnitude for N{ = 3. This suggests 
that including the quarks the perturbative prediction for the potential computed 
through ftqq is also applicable up to a ~ 0.3 and furthermore in full QCD the issue 
of the appropriate scheme is somewhat less important. A direct lattice QCD check 
of these expectations is not possible at present and one has to boldly generalize from 
the Nf — case. In addition, these remarks apply only to the massless theory (mass 
effects were not investigated). Current phenomenological research concentrates on 
application of a velocity dependent potential beyond the static limit for applica- 
tions to top-quark physics [91]. On the one hand, in this application the potential 
is needed for quite short distances, where perturbation theory is intrinsically more 
precise [92], on the other hand, with the velocity entering as a new scale, this rep- 
resents a more difficult multiscale problem. Indeed the renormalization group has 
already been applied to deal with this complication [93]. 

Nevertheless, the lessons learnt in our investigation may be useful in this context as 
well; the type of renormalization group improvement which we found to increase the 
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reliability of perturbation theory has not been applied in [91] so far. 

2.7 Comparison with the bosonic string model 

Although we concentrated on the evaluation of the potential at short and interme- 
diate distances, where the validity of the string picture described in section 2.3 is a 
priori not expected, it is interesting to compare our results with this model. 
If one assumes that the bosonic string picture described is valid already at r = tq, at 
next-to-leading order (NLO) in the effective theory one obtains the parameter free 
predictions 

V(r) = or - + const, (2.90) 
12r 

^0=^ + ^1:5, (2-91) 

where a is fixed by the eq. (2.46), 

ar', = 1.65 - ^. (2.92) 

In Fig. 2.13 one can see that eq. (2.90) (dotted line) is in excellent agreement with 
our results for r > O.Sro. 

As already mentioned in sec. 2.3, in the same region of r, excited potentials do not 
at all follow the expectations from the effective bosonic string theory. This may 
suggest that the agreement with eq. (2.90) is rather accidental. In any case one 
would expect corrections to this formula to be negligible only for much larger r. 
If one includes the higher order correction (NNLO) eq. (2.42) founded in [29] (solid 
line), one can notice that at short distances the agreement is worse, while at our 
largest distances it seems to be better, although our precision here is not good enough 
to make any definitive statement. This fact seems to confirm that the agreement at 
short distances is coincidental. 

The fact that at short distances the NLO prediction is able to describe the data bet- 
ter than the NNLO is not surprising: we are dealing with an asymptotic expansion 
in a region where the expansion parameter ^ is large. Under these circumstances 
it is well known that the introduction of higher order terms in the series do not 
necessarily improve a prediction and can even make it worse. 

Figure 2.14 shows the same comparison for the force; finally, in Fig. 2.15 we 

compared our results for large r with the ones obtained by Liischer and Weisz for 
their largest /? [29]. One observes that the algorithm for the exponential error re- 
duction [64] adopted there is highly necessary to obtain a good precision of the data. 
On the contrary, with our results for r > 1.2ro we would not be able to evaluate the 
correction b in eq. (2.42). 

From this discussion one can conclude that the understanding of the long range 
properties of the static quark potential and the testing of bosonic string model is still 
an important task and requires the employment of advanced numerical techniques 
in order to have a reliable measurement of the observables in that range. 



2.7 Comparison with the bosonic string model 



49 




Figure 2.13: The potential compared with the parameter-free prediction of the 
bosonic string model, fromeq. (2.40) (dashed line) andeq. (2.42) (solid line). The con- 
stant term has been fixed such that the predictions coincide with the non-perturbative 
results at r = 0.9ro. 




Figure 2.14: The force compared with the parameter- free prediction of the bosonic 
string model. The dashed line represents the string prediction at leading order 
eq. (2.91), while the solid line include the correction eq. (2.42) with b = 0.04fm. 
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Figure 2.15: The force compared with the bosonic string model at large distances. 
The points correspond to our computation of the force extrapolated to the continuum, 
while the triangles are taken from [29] for /3 = 6.0. 



Chapter 3 



Scaling properties of RG actions: 
the critical temperature 

Starting from the Wilson formulation of Yang Mills theory on the lattice, eq. (1.5), 
it was observed that for many physical observables it leads to sizeable discretization 
errors. To perform a continuum extrapolation of the physical quantities of interest 
usually requires a big computational effort which is often not practicable. Therefore 
it is highly desirable to use a lattice formulation which reduces the discretization 
errors to have a more reliable estimation of physical quantities. As we have seen, 
the Wilson gauge action shows leading lattice artefacts of order O(a^); by formu- 
lating fermions on the lattice in general one expects corrections starting already at 
0(a). For this reason in recent years most of the efforts were spent to improve the 
fermionic action in view of unquenched simulations. 

For example, by following the Symanzik approach, one can attempt to cancel all 

lattice effects of order a; this development started some time ago [32,94] and the 
coefficients multiplying the various 0(a) counterterms in the improved theory have 
been determined non-perturbatively [95-97]. 

In any case it turns out that also the gauge part plays a very important role, and 
the question which is the most convenient gauge action to adopt has been often 
addressed. 

In particular, the Iwasaki action has been used by the CP-PACS collaboration [37] 
for the most advanced computation of light hadron spectrum. Moreover, the RG 
improved actions have been considered as candidates to be used in the next simula- 
tions on Ginsparg- Wilson/domain wall fermions; interesting features were observed 
for these actions, as the suppression of small instantons and dislocations and a pos- 
sible remedy of the problem of residual chiral symmetry breaking for domain wall 
fermions [98,99]. 

This increasing interest in alternative gauge actions motivates more investigations 
into their properties, starting from the basic ones, like universality and scaling be- 
havior. There are in principle several quantities that one can use to quantify the 
lattice artifacts and to test universality by comparing the results with the plaquette 
action known in the literature. The renormalized observables that we have consid- 
ered are TcTq, where Tc is the deconfining phase transition temperature, and mcro, 
where ma represents the glueball mass and will be discussed in the next chapter. 
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In particular, our study required the computation of the scale r^/a for RG actions, 
which was up to now not present in the literature. 

In this chapter the physical observables will be defined, the details of the lattice cal- 
culation will be explained with particular emphasis on the possible difficulties and 
the results will be discussed. 



3.1 The critical temperature Tc 

In chapter 1 it has been pointed out that Euclidean quantum field theory in the 
functional integral formalism, in particular its lattice regularized version, is formally 
equivalent to a system of classical statistical mechanics in four dimensions. 
Complementary to this, another point of view makes use of Hamiltonian formalism 
and the transfer matrix (see appendix B); from eq. (B.15), eq. (B.18) one obtains 
that for a lattice with Nt points in the time direction and Ng points in the three 
space-directions (we assume for simplicity to be the same for all the space directions) 
^, the partition function of the system is given by 

Z= tr(e-"^*=). (3.1) 

In the limit Ns,Nt — > oo this corresponds to a 4-dimensional quantum field theory 
at physical temperature T — 0. 

For several applications, for example the study of hot hadronic matter or early 
universe, it would be desirable to consider quantum field theory at finite physical 
temperature T also. 

In quantum statistical mechanics the partition function of a finite-temperature sys- 
tem is defined by 

^(^) = tr(e-^«), (3.2) 

where the Boltzmann factor has been set equal to 1, and the thermal expectation 
value of any observable A is given by 

(^)-^tr(e-^«A). (3.3) 

Comparing these expressions with eq. (3.1) and eq. (1.8), one recognizes that a Ng^Nf 
lattice with periodic boundary conditions in the time direction can be regarded as a 
system of finite volume V — {Nsq)^ at finite temperature T, which is related to Nt 
by 

^ = Nta. (3.4) 

The thermodynamic limit corresponds to Ng — > oo; in order to approximate this 
condition in a simulation in a finite spatial volume, one should respect the relation 

Ng > Nt. (3.5) 

^In order to avoid notational confusion, in this chapter we wih use T for the temperature, and 
Nt for the time-extent of the lattice. 
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It was pointed out in [100, 101] that pure Yang-Mills theory undergoes a first order 
phase transition at some temperature T^, below the critical temperature one has 
the confined phase, while the high temperature region is characterized by a gluon 
plasma with freely moving but still interacting gluons. 
The Polyakov loop ^ 

Nt-l 

P{x) = Tr J]^ U{x, 0) (3.6) 

xo=G 

can be interpreted as an order parameter for the deconfining phase transition: in 
the deconfincd phase it assumes a non- vanishing expectation value {P{x)) ^ 0.^ 
In full QCD the low temperature phase corresponds to the usual hadronic QCD 
vacuum, and the high temperature phase is a quark-gluon plasma where the char- 
acteristic low energy features, like confinement and spontaneous chiral symmetry 
breaking arc lost, and the short distance behavior of matter is governed by the 
asymptotic freedom of QCD (see [102] for a review). 

The first non-perturbative lattice determinations of Tc can be found in [103] ; the 
general strategy consists of keeping the relation in eq. (3.4) fixed while varying the 
gauge coupling (3 and therefore implicitly the lattice spacing a. In this way one 
moves the lattice system through the phase transition obtaining 

^ = NM/3c). (3.7) 

There are several methods for determining the critical couphng Pc', it is desirable to 
refer to a definition based on a quantity with a definite finite size scaling behavior, 
for instance the Polyakov loop susceptibility, defined as 

X-V{{\P'\)-{\P\n V^Nl (3.8) 

In the thermodynamic limit, the susceptibility has a delta-function singularity cor- 
responding to a first order phase transition. On a finite lattice the singularity is 
rounded off and one observes a peak for Pc{y)- The critical coupling for each value 
of Nt can be extrapolated to the thermodynamic limit using the finite size scaling 
law for a first order phase transition, 

Pc{Nt,Ns)=Pc{Nt,oo)-hij^j , (3.9) 

where h is considered to be a universal quantity independent of [104] . 

In addition to its intrinsic importance as a fundamental non-perturbative pre- 
diction, Tc provides also a useful quantity to study the lattice artifacts for different 

^Notc that here "Tr" indicates the trace on the gauge group SU(A'^), while in eq. (3.1) "tr" is 
the trace in the Hilbert space. 

^Strictly speaking, only the Polyakov loop correlator {P{x)P^{0)) has a physical meaning and 

can be used as an order parameter. Assuming cluster decomposition at large distances one obtains 
{P{x)P^x)) ~|o|^oo l(-P)l^) ^-iid O'^s can notice that the phase of {P) is not a physically measurable 
quantity. 
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Figure 3.1: Tc/^/a as function of 1/-/V| for different actions, from [50]. By "RG 
improved" is meant here the Iwasaki action. 



Nt 


Wilson [104] 


Iwasaki [38] 


DBW2 [106] 


3 




2.1551(12) 


0.75696(98) 


4 


5.6925(2) 


2.2879(11) 


0.82430(95) 


6 


5.8941(5) 


2.5206(24) 


0.9636(25) 


8 


6.0624(12) 


2.7124(34) 




12 


6.3380(23) 







Table 3.1: The critical coupling from the literature for Wilson, Iwasaki and DBW2 
actions. 

gauge actions and to test universality. 

In order to investigate the scaling properties, one has to measure a dimensionful 
quantity like Tq at zero temperature for (3 = and then build the renormalization 
group invariant combination roTg. For known values of (3c at given for different 
actions one can refer to [38,50,104-106]; for Wilson, Iwasaki and DBW2 actions, 
these values are collected in Table 3.1. 

Instead of using to set the scale, the results that are available in the literature 
are mostly related to the quantity T^j \fd\ they are reported in Table 3.2, where the 
continuum extrapolations are from a reanalysis by Teper [107]. Figure 3.1 collects 
the results for Wilson, Iwasaki, DBW2, Symanzik tree level and FP action, and was 
taken from [50] which reports the latest evaluation of T^j \fo (FP action). 

The first observation from Fig. 3.1 is that from this specific quantity it is not 
possible to arrive at precise conclusions about the discretization errors for different 
actions. 

Furthermore, one notices that for the Iwasaki and the Wilson action a difference of 
order 2(j in the continuum results is observed. The most drastic explanation for this 
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action 



Sym. tree level [104] 
1 loop tadpole impr. [108] 



Wilson [104] 
Iwasaki [38] 
DBW2 [106] 



FP [50] 



oo 0.630(5) 

oo 0.651(12) 

oo 0.627(12) 

oo 0.634(8) 

oo 0.659(8) 

4 0.624(7) 



Table 3.2: Results for the deconfining temperature in units of the string tension 
from different actions [107]. 

discrepancy could be a violation of universality, but this scenario seems unrealistic; 
a more natural explanation is that the string tension is difficult to determine, as 
already pointed out in chapter 2, and it is preferable to use tq to reliably set the 
scale in pure Yang- Mills theory. 

For the Wilson action, the values of ro/a corresponding to the critical couphngs can 
be easily obtained by the parametrization formula [73] ^ 

ln(a/ro) = -1.6805 - 1.7139(/3 - 6) + 0.8155(/3 - 6)^ - 0.6667(/3 - 6)^ (3.10) 

For the Iwasaki and DBW2 action there was up to now no precise evaluation of ro/a 
and we decided to perform new numerical simulations with this purpose. 

3.2 Evaluation of r^/a for RG actions 

For the evaluation of ro/a we followed the same procedure as for the usual plaquette 
action, as described in 2.4.1. We constructed Wilson loop correlation matrices em- 
ploying smearing techniques and multi-hit methods ^. We will not repeat here the 
description of the procedure, but mention only the difficulties with respect to the 
plaquette case. The simulation parameters are reported in Table 3.3 and Table 3.4. 



3.2.1 Effective potential 

We adopted the same smearing operator as for the Wilson action eq. (2.54), with 
M = 4, a = i and the same condition for the smearing levels eq. (2.57). 
For the Wilson action it is possible to define a positive transfer matrix and hence 
a well-defined Hamiltonian; from eq. (2.8) one notices that the contributions to the 
spectral decomposition of correlation functions are always positive. On the con- 
trary, for the Iwasaki and DBW2 action, due to the violation of physical positivity 

^We decided to use the parametrization formula of [73] instead eq. (2.79), which is valid for 
larger range of /3 but less accurate. 

^Due to the presence of (1 x 2) plaquettes in the RG actions, not every temporal link appearing 
in the Wilson loop can be substituted by the one obtained with the multi-hit integration. We 
applied the multi-hit procedure only on alternating links so that two multi-hit links can not belong 
to the same (1 x 1) or (1 x 2) loops. 
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Figure 3.2: The effective potential as function of t for r « ro (ro ~ 4 — 5) for the 
Wilson plaquette action and for the RG-improved Iwasaki and DBW2 action. 



described in appendix B, one expects also negative contributions, so that the plateau 
in the eflFective potential is reached from below. This fact spoils the possibihty to 
find an unambiguous criterion to define the "optimal" smearing: the ability of the 
smearing operator to suppress the higher order excitations can be tested only for 
t ^ tmin, where tmin is estimated in eq. (B.35). 

We started our analysis by observing the time-dependence of the effective poten- 
tial for the different smearing levels. 
Figure 3.2 shows the effective potential 



evaluated from the diagonal elements of the correlation matrices eq. (2.59) for I — 
m — 2, which corresponds to what was estimated to be the optimal smearing for 
the Wilson action. Notice that the negative contributions for the RG actions are 
strongly evident and, as expected, they disappear starting from a certain t which 
becomes larger as the coefficient ci in the action increases. 

One observes however a quite satisfactory plateau in the effective potential, starting 
at sufficiently large t. We decided to extract the potential from eq. (3.11) without 
applying the variational method, which is mathematically founded on the positivity 
of the correlation matrix Cim at a certain t = to; this condition is verified only 
for to ^ tmin, but on the other hand to can not be arbitrarily large because of the 
statistical errors, which increase with t, that make the inversion of Cim impracticable. 
We extracted the potential from eq. (3.11) at t/a — (3 — 4) and we estimated the 




(3.11) 
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systematic error by taking the difference between this value and what one would 
obtain extracting the potential at (t + a) . The systematic errors were linearly added 
to the statistical one; at small r the total uncertainty is dominated by the systematic 
error, while at large distances the situation is in general reversed. 



3.2.2 The force at tree level 

It is interesting to study the force at tree level for different actions, in order to 
investigate how the continuum limit is approached at small coupling ^, where it is 
legitimate to expect that the tree-level approximation is quite accurate. 
Explicitly, the force at tree level is given by [47] 

F,Ur') = ^*-(^) - - (3.12) 

4 Qq n d^k cos{rki/a) — cos{{r — a)ki/a) 

3 a 



, (27r)3 4 j^^3^^ sin'{k,/2) - Ac, ^^'(kj/^)) ' 



Note that for ci = one obtains the results already mentioned for the plaquette 
action eq. (2.50). For the special choice Ci = —1/12, which corresponds to the 
Symanzik tree-level improved action, one finds [47] 



r'V,r(r') = ^| + 0(^)', (3.13) 



if r' = ^/r{r — a); this choice of the action cancels indeed the O(a^) terms at tree 
level. 

In Fig. 3.3 the quantity ^Ftree('"') is plotted as function of {a/r'Y for different 
actions. As expected, in the limit r' ^ oo all curves converge to the classical value 
1^. For the Symanzik tree- level improved action one finds eq. (3.13), that is no 
O(a^) lattice artifacts in the force. One can notice that the RG actions, in particular 
the DBW2, show at tree level large lattice artefacts at small {a/r'). The RG actions 
should reduce the discretization errors at every order in perturbation theory, but it 
turns out that at tree level these actions are "overcorrected" and introduce lattice 
artefacts of the same order or even larger than what is expected with the usual 
plaquette action. This fact should be considered as a problem if one would for 
example evaluate the potential at short distances, as we did for the Wilson action 
in chapter. 2: in particular one should consider the continuum extrapolation with 
great care, because unless a/r is very small, one can not be sure to be in the region 
where the leading discretization errors are quadratic in a. 
For our evaluation of ro/a, the force is defined exactly as in sect. 2.4.2; also in this 
case one can introduce a tree-level improved definition of the force such that 

with no lattice artefacts at tree level. 
From eq. (3.12) one deduces that 

(47rri^)-^ = -[G{r, 0, 0) - G{r - a, 0, 0)]/a, (3.15) 



58 Chapter 3 Scaling properties of RG actions: the critical temperature 



0.13 r~i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

- . Wilson 

■ 

Q 12 — • Symanzik 




0.09 - " - 

- A Iwasaki - 

0.08 ~- o DBW2 

I I I I I I I I I I I I I I I I I 

0.05 0.1 0.15 
{a/rr 



where 

G{r) 



Figure 3.3: The force at tree level for several actions. 



1 r _fk n-=iCos(r,%/a) 

a (27r)3 4(^^^^ sm\k,/2) - Ac, ^^^\k,/2)) ' ^ ' ' 



is the scalar lattice propagator associated to the improved action with coefficient Ci 
in the (1 x 2) plaquette term. Due to the fact that G{f) does not satisfy eq. (C.l) 
because of the presence of higher derivative terms in the action, it is not possible to 
use the coordinate space methods described in appendix C; ri was then computed 
solving the integral eq. (3.16) numerically. 





/5 


ni 


TV 

11 or 


TV 


8 


2.1551 


0,2,4,6 


3 


20000 


12 


2.2879 


0,4,9,13 


4 


4000 


24 


2.5206 


0,12,25,37 


8 


645 


32 


2.7124 


0,18,36,54 


9 


370 



Table 3.3: Simulation parameters for the Iwasaki action for the evaluation of ro- 



3.2.3 Results 

Using a local interpolating formula for the force as explained in sect. 2.4.5, we 
extracted the value of ro/a which are reported in Table 3.5 and Table 3.6. For the 
DBW2 action, besides the three values of Pc known in the literature Table 3.1, we 
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P 


ni 


TV 

ly or 


* meas 


10 


0.75696 


0,2,4,6 


3 


12000 


12 


0.8243 


0,4,9,13 


4 


6000 


16 


0.9636 


0,10,20,30 


8 


800 


24 


1.04 


0,18,36,54 


9 


220 



Table 3.4: Simulation parameters for the DBW2 action for the evaluation of tq. 

decided to evaluate ro/a also for a larger f3 = 1.04, which should roughly correspond 
to P = 6 for the Wilson action and has been used in quenched simulations [99]. Our 
results for the potential and the force at finite lattice spacing are collected in the 
tables A. 5 and A. 6. 





ro/a (r„) 


ro/a (ri) 


2.1551 


2.311(5)(9) 


2.320(6)(9) 


2.2879 


3.026(4) (3) 


3.026(5)(1) 


2.5206 


4.535(6)(4) 


4.511(8)(1) 


2.7124 


6.020(15)(25) 


5.999(15)(19) 



Table 3.5: 

Results for ro/a evaluated at different (3 = (3^ for the Iwasaki action, using the naive 
definition of the force eq. (2.49) or the tree-level improved eq. (3.15). 



(3 


ro/a (r„) 


ro/a (ri) 


0.75696 


2.430(5) (20) 


2.225(4)(11) 


0.8243 


3.129(23)(1) 


3.036(17)(4) 


0.9636 


4.606(13)(17) 


4.556(17)(20) 


1.04 


5.500(29)(7) 


5.452(26)(8) 



Table 3.6: Results for ro/a evaluated at different (3 for the DBW2 action, using the 
naive definition of the force eq. (2.49) or the tree- level improved eq. (3.15). 

The first error contains the statistical uncertainty summed to the systematic one 
due to the interpolation of the force. The second error is the systematic uncertainty 
coming from different choices of t in the effective potential eq. (3.11). One can notice 
that in some cases the last error is the dominating one. 

We can notice that for the DBW2 action the choice of r„ (cq. (2.49)) or ri (eq. (3.15)) 
in the definition of the force leads to results which can be quite different from each 
other, above all at small ro/a; we expected this feature by investigating the force at 
tree level. This ambiguity will make the discussion of the lattice artefacts difficult, 
because the possible conclusions will depend on which definition of the force one has 
used and not on intrinsic properties of the action. 
For the Iwasaki action the two results are not significantly different. 
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Figure 3.4: Parametrization of r^/a for the Iwasaki action, eq. (3.17), using the 
tree-level improved definition of the force. 

3.2.4 Parametrization of r^/a 

Following the strategy of section 2.4.6, one can attempt a phenomenological parametriza- 
tion of ro/a in the range of couplings that we considered. 

For the Iwasaki action the four values of tq (obtained by adopting ri) were fitted in 
the form 

ln(a/ro) = Ci + C2(/3 - 3) + c^{l3 - 3)^ (3.17) 
yielding the numerical results 

ci = -2.1281, C2 = -1.0056, Cg = 0.6041. (3.18) 

in the range 2.1551 < (5 < 2.7124. 

The results and the fit formula are shown in Fig. 3.4; the accuracy is about 0.6% at 
P = 2.1551 and decrease to 0.8% at p = 2.7124. 

For the DBW2 action we have obtained a good accuracy by using a four-parameter 

fit 

ln(a/ro) = di + d2{(3 - 1) + d^iP - i f + d4(3 - l)^ (3.19) 

with 

di = -1.6007, d2 = -2.3179, = -0.8020, d^ = -19.8509, (3.20) 

for the range 0.75696 < /3 < 1.04, where the results always refer to the tree- level 
improved definition of the force. The accuracy is between 0.6% at lower /3 and 0.8% 
at higher /3; it is practically given by the statistical accuracy and the fit is plotted 
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Figure 3.5: Parametrization of ro/a for the DBW2 action using the tree- level im- 
proved definition of the force. 



in Fig. 3.5. 

In [99] the value ro/a = 5.24(3) for the DBW2 action at /? = 1.04 is quoted. Our 
value differs about 3% from that one. 



3.2.5 Scaling of TcVq 

Once ro/a at the given couplings is known, one can finally consider the renormalized 
quantity 

roT„ = i-lH(A), (3.21) 

which has a finite value in the limit a — > (while Nt — > oo). 

The results for Iwasaki and DBW2 actions are given in Table 3.7, together with the 
values obtained using the Wilson action. Also in this case, we show both the results 
obtained with the naive and with the tree-level improved definition of the force. 
The error in T^ro is the quadratic sum of the error for ro/a and the uncertainty in 
(3c, which can be translated to an uncertainty in ro by using the parametrization 
formulas eq. (3.10), eq. (3.17), eq. (3.19), assuming that the coefficients in the fit 
are exact. In our evaluations the error for Pc and the uncertainty in ro are roughly 
of the same order. 

We expect that the leading lattice artefacts are of order a^, such that the continuum 
limit is approached in the following way 



Tero = 7;roU=o + s ■ {aT^f + 0(a7;)^ 



(3.22) 
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Figure 3.6: TcTq for different actions. The x coordinates were slightly shifted for 
clarity. 

The results for Tcro, together with the continuum extrapolation, are shown in 
Fig. 3.6. 

For the Iwasaki action there is no appreciable difference between the results 
obtained with r„ and rj and in both cases the data show better scaling properties in 
comparison to the Wilson action. Furthermore, the value obtained at iV^ = 8 is in full 
agreement with the continuum result evaluated through the Wilson action and hence 
the universality is confirmed; this supports the conclusion that the disagreement 
observed in T^/ \fo is indeed due to the difficulty in evaluating the string tension, 
particularly at small lattice spacings and it is necessary to set the scale through a 
more reliable quantity. 

Also for the DBW2 action the scaling properties are improved, although only using 
Tj instead of r naive, so that it is more difficult to make a statement about the lattice 

artifacts in this case. 

A constrained fit of the form eq. (3.22) including the points with Nt> 6 for Iwasaki 
and Wilson actions yields the continuum result (Fig. 3.7) 

Tcro = 0.7498(50). (3.23) 

At A^t = 6 the Wilson action shows scahng violations for vqTc of about 1.5%, while 
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Figure 3.7: Continuum extrapolation of Tcro for the Iwasaki and Wilson action, 
using the constrained fit eq. (3.22). 



TcVo: 

Nt 


Wilson 


Iwas.(rn) 


Iwas.(ri) 


DBW2(r„) 


DBW2(ri) 


3 




0.7703(50) 


0.7733(53) 


0.8100(90) 


0.7417(90) 


4 


0.7194(22) 


0.7565(23) 


0.7565(23) 


0.7822(65) 


0.7590(56) 


6 


0.7388(30) 


0.7558(33) 


0.7518(31) 


0.7676(68) 


0.7593(76) 


8 


0.7445(33) 


0.7525(60) 


0.7499(55) 






12 


0.7464(44) 











Table 3.7: Results for TcVo- 
they are 0.3% for the Iwasaki action. 

For Nt = 4 the discretization errors for the Wilson action increase to 4%, while for 
the Iwasaki action they remain small (0.6%). 

In Fig. 3.6 we included also the results obtained with the FP action [50], which 

also show a good scaling within 1% even on coarse lattices corresponding to Nt = 3, 2. 
One has however to mention that for those lattices the determination of ro/a contains 
large systematic uncertainties. 

3.3 Scaling of aqq(/^) 

Another interesting observable that can be used to test the scaling violations is the 
dimensionless coupling ctqq obtained from the force. Since in chapter 2 we evaluated 
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the force and the potential very precisely and performed the continuum extrapolation 
for those quantities, we can compare these results with our present determination 
of the force at finite lattice spacings for the Iwasaki and DBW2 actions at different 
distances r. Wc point out that we only determined the on-axis potential and hence 
we can not investigate violations of rotational invariance which would require the 
evaluation of off-axis quantities. 

The coupling ctqq can be obtained from the force by the simple relation 

aqq(A^) = ^F(r)r^ pt^^. (3.24) 

Figure 3.8 shows ftqq in the continuum limit and the results obtained with the Iwasaki 
and DBW2 actions at the largest f] at our disposal. For the Iwasaki action no 
appreciable difference in the results obtained with and ri can be seen, while for 
the DBW2 action the discrepancy becomes large at small distances. 
At large enough distances one obtains a good scaling in the coupling, and one does 
not observe scaling violation within the statistical errors. At small r/ro one sees 
deviations from the continuum limit, as one can observe in Fig. 3.9, where only the 
short distance region is considered. The deviations can be estimated to about 2% 
for the Iwasaki action at r ~ 0.4ro; for the DBW2 action they amount to 4% at 
r ~ 0.4ro if one uses r„ to define the force, and reach even 10% by employing rj. 
This fact shows that the adoption of a tree level improved definition of the force 
does not guarantee success in reducing the lattice artifacts, above all for RG actions, 
which are overcorrected at tree level. 
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Figure 3.8: Qqg at finite lattice spacing for Iwasaki and DBW2 action compared with 
the continuum result. The solid line represents the 3-loop RG perturbative prediction 
of the running coupling as obtained in sect. 2.2; the dashed lines correspond to the 
uncertainty on the A parameter in eq. (2.36). 




Figure 3.9: aqq at finite lattice spacing for Iwasaki and DBW2 action compared 
with the continuum result in the short distance region. 
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Chapter 4 

Scaling properties of RG actions: 
the glueball masses 



Osterwalder and Seller showed in 1978 [63] that in strongly coupled lattice gauge 
theory the lowest eigenstate of the Hamiltonian above the vacuum has a mass m, 
which is usually called mass gap. This result can not be obtained in the framework 
of perturbation theory: on the classical level the gauge field theory does not contain 
any mass term and is scale invariant. Furthermore, the gluon propagator remains 
massless to all orders of the perturbative expansion. This state is hence of purely 
non-perturbative nature (in the infinite- volume limit) and states with higher masses 
are to be expected too. If the masses scale near the continuum limit according to 



it means that they are proportional to the A parameter and hence the correspond- 
ing states survive in the continuum limit and the quantized gauge theory contains 
massive physical states. The existence of these states, called glueballs, has been 
predicted in [109]. 

Being of non-perturbative nature, Monte Carlo simulations on the lattice provide 
a powerful tool to evaluate the glueball masses; they were indeed among the first 
quantities that have been computed on the lattice [110]. 

Since the theoretical discovery of the glueballs, which constitutes one of the most 
fascinating predictions of the pure gauge sector of QCD, also the experimental search 
started. There is indeed an evidence for the existence of exotic glueballs or hybrid 
particles consisting of quarks with gluonic excitation. For a review on the light 
meson spectroscopy see [111]. The difficult task here is to distinguish the glueball 
states from the background of mesons, which have the same quantum numbers. For 
this reason the exotic glueballs, e.g 0"*"", 1"+, are particularly interesting because 
due to their quantum numbers they can not mix with conventional meson states. 
Lattice QCD investigations [112-115], were addressed to study the effects of dy- 
namical quarks and glueball-meson mixing on the glueball spectrum from lattice 



Apart from the physical relevance, the mass of the lightest (O"*""*") glueball is par- 
ticularly interesting since several calculations with the Wilson action [40-43] showed 




(4.1) 



QCD. 
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large lattice artefacts of about 40% at coarse lattice spacings a ~ O.lSfm and still 
20% at a ~ O.lOfm. This quantity could hence in principle provide another stringent 
test on the scaling behavior of alternative gauge actions. It is desirable to check the 
discretization errors on several physical observables, since they can be quite different 
depending on which quantity we are considering. 

As already pointed out since the first works, the calculation of glueball masses 
presents a lot of technical difficulties; due to the fact that these are quite large 
{ma > 1.6GeV), the signal in the correlation functions of the gluonic excitations 
decay fast and disappears in the noise. 

Several smearing techniques were developed in order to suppress excited states and 
to be able to extract the masses at moderately large value of t. It turned out that 
these techniques are more efficient in the region of small lattice spacings, where one 
expects for this reason the determination of the glueball masses to be easier. 

However, as already pointed out, one has to remember that in order to neglect 
finite size effects one has to consider a physical lattice extent L > 3ro, so that the 
simulations at small lattice spacing arc very expensive and it is very difficult to reach 
a large statistics which is usually required for those calculations. 
A possible solution to this problem can be to use anisotropic lattices [116], where 
the resolution in the time direction is finer that in the space directions <S a^, and 
hence one can follow the signal over a larger range of time slices. Morningstar and 
Peardon [117,118] applied this method to evaluate the glueball spectrum below 4GeV 
in the pure SU(3) gauge theory and improved the previous determinations of the 
glueball masses remarkably. In this chapter we will present the calculation of amo++ 
and am2++ for RG actions in the range of lattice spacings 0.11 fm ^ a ^ 0.17fm, 
and then make use of the knowledge of rg/a to build the renormalized quantities 
romQ++ and rom2++ . 

Our intention here is not to have a more precise determination of the glueball masses 
than what can be found in the literature and indeed our total errors on the measure- 
ments are quite sizeable; we will focus on the investigation of discretization effects 
and on possible further difficulties that can arise in the extraction of the masses from 
correlation functions. 



Physical states in the vacuum sector of the Hilbert space of lattice gauge theory can 
be obtained by applying gauge invariant operators to the vacuum. As particular 
candidates one can take space-like Wilson loops lying in the xo = timeslice in the 
fundamental representation eq. (2.2), considered as multiplicative operators on wave 
functions. The simplest choice is the plaquette eq. (1.4) 



4.1 Glueball states 




(4.2) 
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where x = {xi,X2, x^), i, j E {1, 2, 3}. 

Applied to the vacuum the operator yields a state 

l*..-(^)> = {Wl/^'ix) - (0|W,f ^(f)|0)} |0), (4.3) 

where the projection onto the vacuum itself has been subtracted. 
Eigenstates of the momentum operator are formed by Fourier transformation 

|*,,(^)) = L-^/2X^e-^-^^|^,,(f)). (4.4) 

X 

If one is interested in the determination of masses it is sufficient to consider zero 
momentum states, denoted by 

I*,,.) = |*,,.(0)). (4.5) 

We assume that the continuum limit of SU(3) lattice gauge theory exists and their 
rotational invariancc is restored: the spin of a state is then characterized by the 
unitary irreducible representations of SO (3), for bosonic states. 
On a cubic lattice the rotation symmetry is broken down to the cubic group O; the 
eigenstates of the Hamiltonian have hence to be classified according to the unitary 
irreducible representations of O. 

For the first application of the representation theory of the cubic group to lattice 
gauge theory one can refer to the paper by Berg and Billoire [119]. 
The cubic group contains 24 elements, corresponding to the permutations of the 
four space diagonals of a cube. There are five irreducible representations, denoted 
by Ai, A2, E, Ti and T2, which have dimensions 1,1,2,3 and 3. Since the cubic group 
is a subgroup of SO (3), any representation Dj of SO (3) for a state of spin J will 
induce a representation on the group O, the so-called subduced representation Dj. 
It will in general no longer be irreducible and can be decomposed into the irreducible 
representations of O. Up to J = 4 one finds 



= A (4.6) 

= Ti (4.7) 

= E®T2 (4.8) 

Df = AaeTieTa (4.9) 

d2 = Ai®E®Ti®T2. (4.10) 



For example, a spin 2 particle is described in the continuum by a quintuplet of de- 
generate states; on the lattice the quintuplet is split into a doublet E and a triplet 
T2. At finite lattice spacing one expects a mass splitting between the two represen- 
tations, so that ruE ^ 171x2- By approaching the continuum limit the ratio mE/mT2 
will converge to one in order to restore the full Euclidean rotational symmetry. For 
the Wilson action we expect that the leading discretization corrections in mass ratios 
are of order ; for other lattice formulations this is no more obvious and in fact with 
our data we will not be able to give conclusions on the form of the lattice artefacts 
in the specific range of lattice spacings that we have considered. 
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Given a lattice operator which transforms according to some irreducible representa- 
tion R of the cubic group, by applying it to the vacuum one obtains a state 

i*ii) = E^«i*«)' (4.11) 

a 

where \^a) are eigenstates of the Hamiltonian. In the continuum limit any \^a) 
belongs to some spin J multiplet, so that \^ r) contains various spins J. A spin J 
can occur in this superposition only if R is contained in . 

At this point the important assumption is that the lowest spin contained in R will 
belong to the lowest mass and will therefore dominate the correlation functions. The 
representations of the cubic group, their dimensions and the lowest spin content are 
summarized in Table 4.1. 

In addition to the symmetry associated with the cubic group, there are two discrete 
symmetries: the total space reflection and the charge conjugation. The eigenvalues 
of the space reflections correspond to the parity P = ±1; the cubic group combined 
with space reflection forms the group O/^ = O x Z2, which contains 48 elements. 
The charge conjugation corresponds for the Wilson loops to the complex conjugation, 
and its eigenvalues are C = ±1. The states belonging to an irreducible representation 
of the lattice symmetry are then labeled by 

I*) = |J^^). (4.12) 

4.2 Irreducible representations of the cubic group 
on Wilson loops 

Spatial Wilson loops, being gauge invariant operators acting on the vacuum, create 
physical states like eq. (4.5). Here we will consider the space-like Wilson loops up to 
length 8, for which all irreducible representations of Oh have been constructed [119]. 
The loops are shown in Fig. 4.1. 

The first simple rule that one can easily obtain is that real parts of the Wilson 
loops have C-parity C — imaginary parts have C — —1. Every Wilson loop 
corresponding to a path of length L can be represented by an L-tuple 

L 

{h,...jL), With 5]/i = 0, (4.13) 

i=l 

where fi are unit vectors corresponding to the space-like coordinates in positive and 
negative directions (±1, ±2, ±3). The equivalence class corresponding to (/i, fi) 
will be denoted by [(/i, Jl)] and will contain L-tuples which are equivalent up to 
cyclic permutations. 

Under C-parity one has the transformation 

^[(/i,-,/l)] = [(-/l,...,-/i)], (4.14) 

so that the combinations 



[(A, /l)]± = [(A, A)] ± [(-A, -A)] (4.15) 
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R 


Dimension 


Lowest spin J 


A, 


1 





A2 


1 


3 


E 


2 


2 


Ti 


3 


1 


T2 


3 


2 



Table 4.1: Irreducible representations of the cubic group O. 
are even and odd under C-parity respectively. 

On operators of fixed shape we generate a representation M. of Oh by means of 

MA(fi,...,fL)]±^[{Mji,...,MjL)]±, ygeOh, (4.16) 

where M.g is the matrix corresponding to the group element g in the fundamental 
representation. Excluding the C-parity, the dimension d of the generated represen- 
tation will be less or equal than 48, which is the order of the group Oh and which 
will correspond to the number of different spatial orientations of the given shape. In 
the second column of Table 4.2 the dimensions d of the generated representations 
on every loop shape are listed. 

The irreducible contents of the representation M. can be determined by means of 
the so-called character relation. The characters are a set of quantities which are the 
same for all equivalent representations. For finite groups they uniquely determine 
the representations up to equivalence, providing in particular a complete specifica- 
tion of the irreducible representations contained in a given representation. In general 
we will have 

M ^mr^e+nar^e..., (4.17) 

where are the irreducible representations listed in Table 4.1 and Up is the multi- 
plicity of F^ in 
For fixed C-parity one finds 

np^^Y^ncxiWiC), (4.18) 

c 

where the sum goes over all classes of conjugate elements, ric is the number of 
elements in the class C and x(C), x^(C) are the characters of F^ respectively; 
for matrix representations the character is given by the trace of the corresponding 
matrix and is the same for all elements in a given conjugacy class. The results for 
the C-parity plus representations are given in Table 4.2. 

4.3 Wave functions of glueball operators 

In order to perform a Monte Carlo calculation of the glueball masses, for each ir- 
reducible representation one has to calculate an orthonormal basis explicitly. This 
procedure involves the character projection operator, which is defined in general by 

^ Teg 
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#1 



#2 
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#4 
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Figure 4.1: Wilson loops un to length 8. 
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where p labels the irreducible representations of dimension dp of a finite group of 
coordinate transformations Q of order g, x^(T) being the character of T G ^ in 
and P{T) the unitary operator in the Hilbert space of the coordinate transforma- 
tions T eQ. 

projects out of a function (f) of the Hilbert space the sum of all the parts transform- 
ing under T^. Having chosen a function such that V^(f) does not vanish identically, 
one constructs P{T){V^4') for each T & Q. Each of these are linear combinations of 
the dp basis functions of F^. Prom these functions one extracts dp linearly indepen- 
dent functions and applies the Gram-Schmidt orthogonalization to obtain a set of 
orthonormal functions of the irreducible representation. 

What we practically did is to construct from a loop prototype for each shape by 
respecting the following rule: one builds the path first going in 1-, then on 2- and 
finally, if necessary, in 3-direction. Prom this reference orientation all others can 
be generated by applying the group elements of Oh in a given order, taking care 
not to generate orientations that are equivalent up to translations. Once all the 
independent orientations were created, we formed the orthonormal basis systems 
for each irreducible representation by taking the related coefficients from the litera- 
ture [119,120]. 

As example we can consider the single plaquette operator (loop shape #1), for which 
the orthogonal wave functions are listed in Table 4.3, with the following convention. 
The numbers in each row denote the contribution of the specific orientation to the 
wave function in question; the three positive orientations of the single plaquette can 
be labeled as O12, O13 and O23. The wave function associated to the irreducible 
representation A'^'^ will be constructed by forming 




O12 + Oi3 + O23, 



(4.20) 



while the E'^'^ functions are 



O 



,E++ 




(4.21) 
(4.22) 



1 



23 



Of 



,E++ 



Oi3 — 023- 
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shape 


d 


A++ 


A++ 


Hi 




rp + + 


^1 




77"+- 
11/ 


^1 




#1 
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1 
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#2 
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2 























#3 


12 
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#8 


48 
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#17 
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1 





1 
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#18 


24 
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1 
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1 
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#19 


6 


1 





1 
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#20 


24 


1 





1 
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1 
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#21 


12 


1 
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1 
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#22 


48 
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2 
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Table 4.2: Irreducible contents of the representations of the symmetry gronp of 
the cube on Wilson loops up to length 8, for C = +1. d is the dimension of the 
representation of Oh on the different loop shapes. 

4.4 Simulation details 

For the numerical simulation, we followed essentially the method proposed by [50]. 
We performed a test simulation using the standard Wilson action for /3 = 6.0 on a 
16^ lattice, for which rcsTilts can be found in the literature. 

We decided to concentrate on mo++ and m2++ by measuring the masses in the repre- 
sentations Ai~^, E^^ and Tg*"^. The mo++ is particularly interesting for our purpose 
to investigate the lattice artefacts on RG actions since these turn out to be quite 
sizeable for the Wilson action. 

For our test simulation we measured all 22 loop shapes up to length 8 and formed 
the wave functions corresponding to Al^ , E^^ and as explained in the previous 
sections. Our timeslice observable at time t is 

r-3/2 f^ri 

S^it) = ^Y.Y.'^^^n{^^^^ ^ = l,-,22, (4.23) 

X i=l 

where the coefficient are taken from the literature [119,120], the sum over i 
indicates the sum over all dn orientations of a given shape n and X is a suitable 
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loop shape #1 
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1 
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-1 
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Table 4.3: Orthogonal wave functions of the irreducible operators that can be built 
from the plaquette. 

normalization constant. 

As for the calculation of the potential, a smearing procedure is necessary in or- 
der to improve the overlap with the ground state. We applied the smearing op- 
erator eq. (2.54) to every space-like link; for j3 = 6.0 we found that M = 4, 
ni = (5, 10, 15, 20) was a reasonably good choice for the smearing levels. 
Then we build the correlation matrices ^ 

Ciit) = {S^{t)SPm^ = {S^{t)S['m - {S§{t)){S^{0)), (4.24) 

for R = T^^. The indices k, I assume (22 x M) values, where M is the 

number of smearing levels. 
We expect that 

{S^{t)SP{Q)), = |e-^*|*f) = |*„)(*«|*f)e-''(^)*, (4.25) 

a 

where {"^a) are eigcnstatcs of the Hamiltonian. At large t the lowest mass mQ{R) 
dominates; it belongs to a glucball state which in the continuum limit will have the 
lowest spin contained in the representation R. 

Starting from our large basis, we analyzed the signal/noise ratio of the different 
operators in order to eliminate those which introduce large noise and to reduce the 
correlation matrix to a set of well measurable operators. It is important not to 
include noisy operators in the analysis performed through the variational method, 
because they could introduce numerical instabilities in the generalized eigenvalue 
problem. 

We considered the effective masses evaluated from the diagonal elements of the 
correlation matrix eq. (4.24) 

^Notice that the vacuum subtraction is required only in the A^~^ channel, since it has the same 
quantum numbers as the vacuum. 
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and took the relative error as quality of the corresponding operator. 

We found general agreement in the classification of the "bad" operators with [120]. 

Our final choice for the operators to measure was then 



#2, #5, #7, #8, #10, #12, #14 
#2, #5, #7, #10, #12, #14 
#7, #12, #14, 



(4.27) 
(4.28) 
(4.29) 



where for and T^'^ channels we took respectively 2 and 3 orthogonal projec- 
tions for each shape. For each operator we then considered all M smearing levels. 
One has to keep in mind that our choice contains a certain degree of arbitrariness, 
since the behavior of the different operators may depend on the considered lattice 
spacing. We did not perform a systematic study and we assumed that this choice 
is reasonable for each value of the couphng and also for the Iwasaki and DBW2 
actions. 

The simulation parameters for our new Monte Carlo simulations are listed in Ta- 
ble 4.9. Measurements were taken after a number of sweeps between three and five. 



4.5 Analysis details 

For the extraction of glueball masses from the correlation matrices eq. (4.24) varia- 
tional techniques are necessary. Due to the violation of physical positivity for the RG 
actions (appendix B), the variational method is mathematically not well founded, 
at least not at very small to, where one would like to apply it. The statistical errors 
are indeed drastically increasing with t and already at t = 4a the signal in the cor- 
relation function eq. (4.24) is lost in noise. 

We decided to adopt the following procedure: first we solved the generalized eigen- 
value problem 

C{h)Vaitl, to) = XaC{to)Va{tuto), (4.30) 

with to — 0, ti — a. Then we projected the correlation matrices to the space of 
eigenvectors corresponding to the N eigenvalues which satisfy the condition ^ 

A^>e, (3 = 1,...,N, (4.31) 

where e is an adjustable (small) parameter. So the reduced matrix is obtained by ^ 

C^{t) = {v,{ti,to),C{t)vj{ti,to)), t,3 = 1,...,N. (4.32) 

By choosing e appropriately in eq. (4.31) one hopes to get rid of most of the unphys- 

ical modes caused by negative and very small eigenvalues. 

Then one can apply the variational method to the reduced matrix 

C^(i)w^(i, to) = A^(i, to)C''{to)wp{t, to). (4.33) 

^in [50] the direct eigenvalues and eigenvectors of C (to) are considered, instead of the generalized 
one. We have tried both possibilities and found consistent results. 
^The index R for the representation is now omitted. 
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The effective giuebaii mass can be read off directiy from tlie lowest eigenvalue 



or, alternatively one can project again C to the subspace corresponding to the 
lowest eigenvalue 



We tested that these two different evaluation yield results which are compatible 
within the statistical errors. 

We chose e in cq. (4.31) so that the reduced matrix had dimension between 2 and 5. 
We applied the same procedure also for to = a, ti = 2a in eq. (4.30), but due to 
the fact that the statistical fluctuations are already quite large one has to start from 
the beginning from a reduced number of operators in order to be able to solve the 
generalized eigenvalue problem. 

We decided to extract the glueball masses by taking our effective mass at t = 3a for 
the small (3 regime and t = for the large value of the coupling at our disposal. For 
the RG actions we performed numerical simulations up to L = 16, with a minimum 
lattice spacing a ~ 0.1 fm; in this sense we are quite far away from the continuum 
limit and we will not be able to perform a continuum extrapolation of the results, 
but on the other hand this is the regime of lattice spacings that have been used for 
simulations with dynamical fermions and hence it is desirable to obtain informations 
about the discretization errors in this region. 

The figures 4.2,4.3 show the effective masses in the Af'^ channel computed by using 
eq. (4.34) (filled squares) and eq. (4.36) (empty squares) for the Iwasaki and the 
DBW2 action, with to = 0. For the largest value of (3, which are more precise, we 
show also the results obtained with to = a (filled triangles). As for the potential, 
one notices also in this case the presence of negative contributions in the correlation 
functions due to unphysical states; the situation here is somehow more drastic since 
one has to discard for this reason the small t region but on the other hand the errors 
increase very rapidly and one is forced to extract the glueball mass ait = 3 — 4a. The 
results can then be affected by systematic errors that can not be easily estimated. In 
the plots we have indicated with the dotted lines the range in which we decided to 
take the mass. The statistical errors were evaluated by using a jackknife procedure. 
Concerning the determination of the m2++ we observed that the signal for the E~^^ 
channel is usually worse than for the T2~^ and the errors on the effective masses are 

^Due to the periodic boundary conditions in the time-direction, the correlation matrix eq. (4.24) 
describes not only the correlation C{t), but also C(T — t), where T is the lattice time-extent. 
Therefore one would expect C{t) ~ ^(e-"** + e"™'^"*)). If one assumes that T is much larger 
than the value of t at which we are able to evaluate the mass, one can neglect the contribution 
from {T — t). For our largest lattice spacing (T = 10a) we included this corrections and check that 
they are irrelevant for the measurement of the glueball masses. 




(4.34) 




(4.35) 




(4.36) 
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very large already at t = 3a. For this reason we decided to use mrp++ as estimate of 
1712++ at finite lattice spacing. 

The figures 4.4,4.5 show the effective masses for the T^'^ channel. The determination 
of the masses is more problematic that the Af'^ case; in particular for the largest 
coupling one can observe that at t = 4a the errors are too large to have a significant 
measurement of the mass. Wc decided hence to extract the mass at t = 3a, taking 
care that this value is compatible with the one at t = 4a within the statistical errors. 



4.6 Results 

In Table 4.5, Table 4.6 the results for the masses in lattice units are reported. For 
the smallest f3 for the Iwasaki action, it was not possible to obtain a reliable evalu- 
ation of the mass in the T.^^ channel. 

The renormahzed quantities roma are reported in Table 4.7, Table 4.8; due to the 
fact that the errors on amc are quite large (between 8% and 9% for Ai~^ and be- 
tween 10% and 20% for T2^^), the difference which arises by choosing rj or r„ is in 
this case much smaller than the total uncertainty; we have reported however both 
results in the table. For (3 = 2.2423 we had no direct measurement of ro/a and we 
used the interpolation formula eq. (3.17) to evaluate it. 

Figure 4.6 shows the results for romo++ as function of (a/ro)^, displaying only the 
results obtained with r„. For the comparison we included the results for FP ac- 
tion [50] and several calculations performed with the Wilson action [40,116,121]. 
The continuum values available in the literature are listed in Table 4.4 and have 
been taken from [50] for the FP action and from [122] for the Wilson action, where 
the results of [40, 116, 121] have been expressed in units of tq. 
The interpretation of our results is not very clear, also due to the large errors. 
In any case, we want to stress that our determination can be seen at least as an 
upper limit for mo++ and m2++. We expect that at the values of t/a at which we 
extracted the masses, the effects of positivity violations have already disappeared; 
this assumption is justified by our study based on perturbation theory and our esti- 
mation of tfriin eq. (B.35). 

For this reason we believe that possible systematic uncertainties on the glueball 
masses could only be due to the presence of excited states and hence could affect 
our measurement only in such a way that the real values of mo++ , m2++ are lower 
with respect to our determination. 

At lattice spacings a ~ 0.15 fm we notice a improvement of the RG actions with 
respect to the Wilson action; comparing with the continuum limit we find no signif- 
icant discrepancy both for DBW2 and Iwasaki action, while for the Wilson action 
one finds 30 — 40% deviation. 

At lattice spacing a ~ 0.1 fm we find on the other hand large lattice artefacts for RG 
actions: the result obtained with the Iwasaki action is practically compatible with 
the one calculated through the plaquette action at the same lattice spacing, while 
for the DBW2 action it is even further away from the continuum limit. 
If one considers our measurement as upper limit, one could conclude that the RG 
improved actions are not able to cure the problem of large lattice artefacts for the 
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Figure 4.2: The effective masses for the A^~^ channel, evaluated with Iwasaki action, 
at different lattice spacings. The filled and empty squares corresponds to respectively 
to eq. (4.34) and cq. (4.36) with to = 0. For /3 = 2.5206, the filled triangles correspond 
to eq. (4.34) with to = a. 
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Figure 4.3: The effective masses for the A^'^ channel, evaluated with DBW2 action, 
at different lattice spacings. The filled and empty squares corresponds to respectively 
to oq. (4.34) and eq. (4.36) with to = 0. For (3 = 0.9636, the filled triangles correspond 
to eq. (4.34) with to = a. 
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Figure 4.4: The effective masses for the T^'^ channel, evaluated with Iwasaki action, 
at different lattice spacings. The filled and empty squares corresponds to respectively 
to eq. (4.34) and eq. (4.36) with to = 0. 
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Figure 4.5: The effective masses for tlic T2 cliannel, evaluated with DBW2 action, 
at different lattice spacings. The filled and empty squares corresponds to respectively 
to eq. (4.34) and eq. (4.36) with to = 0. 

0"'"''" glueball mass. At very small lattice spacings one expects that the dominant 
lattice artefacts arc of order a^; for the Wilson action this is indeed well confirmed 
by the numerical results in the range a < 0.17fm, as one can see in Fig. 4.6. 
For alternative actions there is no reason a priori to observe the same behavior: 
while at small lattice spacings the O(a^) should in any case be the dominant one, 
at larger a it is possible that lattice artefacts are governed by higher orders (O(a^) 
and higher). 

We have indeed already observed deviations from O(a^) behavior in RG actions for 
the quantity T^ro and even already for the force computed at tree level (chapter 3). 
With our results it is not possible to investigate how the continuum limit is ap- 
proached, because one should evaluate the glueball masses at smaller lattice spacing 
and this was beyond the purpose of this work. The evident problem shown in our 
plot is that a continuum extrapolation of the form 



romo++ = romo++\a=o + s x 




(4.37) 



of the data obtained with the RG action in the considered range of lattice spacings 
would lead to results which are rather different from the expected continuum limit. 
One has however to remember that the extraction of the glueball masses is for RG 
actions particularly difficult and the systematic errors could be underestimated. 

In Fig. 4.7 we report our results for rom2++; for this particular observables the 
calculation performed with the Wilson action do not show significant lattice arte- 
facts. 

At our smallest lattice spacing we do not observe a deviation from the results ob- 
tained with the Wilson action; one has however to notice that our errors are too 
large to make any conclusive statement. 
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Collab. 


romo++ 


rom2++ 


M & P [118] 


4.21(11)(4) 


5.85(2)(6) 


GFU [121] 


4.33(10) 


6.04(18) 


Teper [107] 


4.35(11) 


6.18(21) 


UKQCD [40] 


4.05(16) 


5.84(18) 


FP [50] 


4.12(21) 


[5.96(24)] 



Table 4.4: Continuum extrapolations of the two lowest glueball masses in units 
of To- For the FP action, the 2++ value is not extrapolated to the continuum but 
denotes the mass obtained at a lattice spacing a = 0.10 /to. 



Iwasaki action 



p 




amrp++ 


2.2423 


1.11(10) 




2.2879 


1.20(7) 


1.50(32) 


2.5206 


0.72(6) 


1.17(10) 



Table 4.5: Glueball masses in lattice units, Iwasaki action. 



For our largest lattice spacing the results with the Wilson action are not available 
(there are only know results with anisotropic lattices [118]) and it is indeed difficult 
to have a rehable estimation of the mass also in our case: we decided to show how- 
ever our results, even with the very large error bars. 

From our computation one can not deduce if for rom2++ the RG actions show sig- 
nificative discretization errors and further investigations are needed to clarify the 
issue. 



DBW2 action 



p 






0.8243 
0.9636 


1.27(10) 
0.62(7) 


1.70(39) 
1.22(15) 



Table 4.6: Glueball masses in lattice units, DBW2 action. 
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Iwasalvi action. ?•„ 



p 


rom^++ 


romrp++ 


2.2423 

2.2879 
2.5206 


3.08(28) 

3.63(21) 
3.26(27) 


4.54(97) 
5.31(45) 


Iwasaki action, rj 




rom^++ 


rom^++ 


2.2423 
2.2879 
2.5206 


3.07(28) 
3.63(21) 
3.25(27) 


4.54(97) 
5.28(45) 



Table 4.7: Results for romo for the channels Af~^ and T2"'', using the Iwasaki 
action. 



DBW2 action, r, 







romrp++ 


0.8243 
0.9636 


3.97(31) 
2.86(32) 


5.3(1.2) 
5.62(69) 


DBW2 action, n 




rom^++ 




0.8243 
0.9636 


3.86(30) 
2.82(32) 


5.2(1.2) 
5.56(68) 



Table 4.8: Results for romo for the channels Af~^ and T2"'', using the DBW2 action. 



Iwasaki action 





L 


ni 




2.2423 


10 


2,4,6,8 


12000 


2.2879 


12 


2,4,6,8 


16000 


2.5206 


16 


3,6,9,12 


8000 


DBW2 action 


P 


L 






0.8342 


12 


2,4,6,8 


8000 


0.9636 


16 


3,6,9,12 


2500 



Table 4.9: Simulation parameters for the evaluation of the glueball masses for the 
Iwasaki and DBW2 actions. 
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Figure 4.6: The 0++ glueball mass normalized with ro as function of (a/ro)^ for 
different actions. 
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4.7 Open questions 

As already mentioned, RG action have been indicated by some Collaborations [37, 
98, 99] as best candidates to be used in the simulations with dynamical fermions. 
In view of these next simulations, it is important to have under control the proper- 
ties of alternative gauge actions by testing their features, the universality and the 
discretization errors for several physical quantities. 

This motivated our investigations and requires further studies on this subject. 
An important point to stress is that in full lattice QCD it can happen that the fea- 
tures of different gauge actions are modified, also depending on which formulation 
of fermions on the lattice is adopted. 

For example, the JLQCD and CP-PACS collaborations have used a 0(a) formula- 
tion of the fermionic action, the so-called Sheikholeslami-Wohlert (clover) action [32], 
which is based on the Symanzik approach (see chapter 1). In order to achieve the 
improved action one has to add a countertcrm such that the order a term in the 
action of the effective continuum theory is canceled. For the improved action one 
thus obtains 

S = Swiison + «^ X] Csy,ilj{x)^a^^Fi^^{x)ilj{x), (4.38) 

X 

where Swuson corresponds to the Wilson fermionic action ^ and F^,^ is a lattice rep- 
resentation of the gluon field tensor. 

The coefficient Csn, in cq. (4.38) has to be tuned appropriately. For example, for 
simulations with 2 flavors of dynamical fermions, the CP-PACS collaboration used 
the mean flcld 0(a) improved definition for the fermionic action together with the 
Iwasaki gauge action [37], while the JLQCD collaboration used non-perturbative 
0(a) improved Wilson fermion together with the plaquette gauge action [123]. 
It was observed [124] that for Nf — 3 simulations with 0(a)-improved Wilson 
fermions and plaquette gauge action one finds an unexpected first order phase tran- 
sition at zero-temperature. 

The investigations indicate that this phase transition is a lattice artefact restricted 
to strong couphng regions (/3 < 5.0). This phenomenon spoils the possibihty to 
perform a continuum extrapolation of physical quantities unless one considers only 
couplings in the region f3 > 5.0. 

The interesting fact is that this transition disappears if one adopts alternative gauge 
actions, like Iwasaki or tadpole-improved Symanzik action, or alternatively if one 
uses the unimproved Wilson fermionic action (cgw — 0) with the plaquette gauge 
action. 

A plausible scenario is that the clover term added in eq. (4.38) produces a term in 
the adjoint representation of the gauge fields with a positive coupling. 
It is indeed known that pure SU(3) lattice gauge theory containing fundamental and 
adjoint representations of the gauge fields undergoes a first order phase transition 

^Fermions on the lattice are described by Grassmann variables tp, ip defined on the lattice 
sites. There are several possibilities to define the lattice action associated to fermions, leading to 
different properties. The so-called Wilson action on the lattice can be written in the form Swuson = 
'4'{x){D + mo)V'(a^), where mo is the bare quark mass and D = 5{7^(V* + V^) — aV* V^} 
is the lattice Dirac operator. 
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for positive values of the adjoint coupling [125-127]. A simple way to construct a 
gauge action with an adjoint coupling would be to generalize the Wilson plaquette 
action in the following way: 

K^^'i^)] + /^^ E E {i - l\K^'i^)\"] ' (4-39) 

where W^^^ is the plaquette defined in eq. (1.4). 

The resulting line of the first order phase transition with the location of the end- 
point for the gauge group SU(3) is shown in Fig. 4.8. The Wilson plaquette action 
corresponds to the line 0^ = 0: it is located below the end point of the phase tran- 
sition and hence it is safe to take the continuum limit along this line. If one would 
approach the continuum limit by crossing the transition line above this point, there 
would be discontinuities in the expectation values of Wilson loops and in glueball 
masses. In particular, at the point of intersection with the transition line the 0"'""'" 
glueball mass is expected to be zero. 

The fact that for the Wilson action the lattice artefacts for mo++ are large and the 
glueball mass in lattice units becomes small could then also be interpreted as the 
"influence'' of the endpoint of the first order phase transition on the Pa — ^ line 
corresponding to the Wilson action. 

In order to test the validity of this scenario, one should add an adjoint term in the 
action and check how does it change the behavior of the 0"*"+ glueball masses. 

An important test is also the computation of scalar glueball mass in full QCD 
^; recent results [128] are shown in Fig. 4.9. The SESAM [114] results have been 
obtained with Wilson gauge and fermion action, while UKQCD [115,129] adopted 
0(a) Wilson improved fermions. The larger lattice spacing results show a significant 
reduction in the lightest scalar mass, as shown in Fig. 4.9. 

It is very important to understand if it implies a lower scalar mass also in the con- 
tinuum limit or if an enhanced 0{a?) correction might be present. 

A useful study would be to repeat these investigations with alternative gauge 
actions, starting for example by locating the position of the end point of the first 
order phase transition; it is interesting to know if it lies closer or further away from 
the (3a = ^ line with respect to the plaquette action. 

From our calculation of the glueball masses with Iwasaki and DBW2 action one 
can not deduce any possible scenario and a systematic investigation which includes 
adjoint terms would be needed. 

Discussing these open issues will be helpful to understand properties of full QCD 
on the lattice and choose an optimal combination of gauge and fermionic actions to 
perform reliable unquenched simulations. 



^By inserting dynamical fermions in the theory, the mass eigenstates with 0^"*" quantum num- 
bers are not distinguished as "gluebaUs" or as "quark-antiquark" states; in general mixing will 
occur. 
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Figure 4.8: The phase diagram of SU(3) gauge theories in the fundamental- adjoint 
couphng space [130]. 




Figure 4.9: Plot from [128]for the value of mass of the 0++ glueball state from 
quenched data [40-42, 118] and for the lightest scalar meson with Nf = 2 flavours of 
dynamical quarks [114, 115, 129]. 
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The main object of the first part of this work was the potential between a static 
quark and antiquark in the pure SU(3) gauge theory. In particular, we concentrated 
on the question whether at short distances large non perturbative terms are present. 
To test this possibility we used both the potential and the force, which define two 
different renormalization schemes. The first step was to obtain a parameter-free 
perturbative prediction for these quantities from the Renormalization Group equa- 
tion at two and tree loops; this was possible due to the fact that the A parameter 
is known from lattice studies [15] and the coefficients of the /3-function have been 
calculated up to tree loops [20,53-57]. 

Then we performed a numerical Monte Carlo simulation to evaluate potential and 
force non-perturbatively. In the pure gauge theory it was possible to compute these 
quantities over a large range of distances with high precision by considering two sets 
of lattices and matching them (in the continuum limit) in a region of overlap. By 
setting the length scales tq and Vc we built renormalization group invariant quanti- 
ties and we performed a continuum extrapolation of our results. At this stage we 
investigated the scaling violations and we were able to confirm Symanzik's picture 
of lattice artifacts. By comparing these results with the perturbative prediction we 
did not observe the presence of large non-perturbative terms. Furthermore, we re- 
marked that different renormalization schemes have different convergence properties 
and the scheme defined through the force is in this sense more appropriate: here 
perturbation theory works remarkably well up to ~ 0.3, while for the potential 
only up to ay ~ 0.15. 

This result confirms what can already be argued from the apparent "convergence" by 
comparing the two- and three-loops expansions for the running coupling in the two 
schemes; the fact that the perturbative coupling defined through the potential fails 
in describing non-perturbative data is due the restricted domain of applicability of 
perturbation theory for this scheme rather than a signal of non-expected terms. 
Although the situation in full QCD is much more complicated due to presence of 
different mass scales, this study may be a guideline for future investigations and 
applications of perturbation theory. 

We then compared our results with the bosonic string model, which is expected 
to describe the long-distance properties [26] of the theory. We observed a surprising 
good agreement with the next-to-leading order prediction [27, 28] of the effective 
theory already at 0.4fm, where in principle the validity of the string picture is ques- 
tionable. 

By including a higher order correction estimated in another work [29] we observed 
that at large distances the agreement with the non-perturbative data is improved. 
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while at intermediate distances it is deteriorated. This fact is not surprising since 
we are deahng with asymptotic expansions and it could indicate that the agreement 
at intermediate distances is rather accidental. However we want to remark that we 
evaluated the potential only up to ~ 0.8fm, but a deep understanding of the string 
nature of the strong interactions requires the investigation at larger distances, where 
the statistical errors are usually very large. With the help of the numerical tech- 
niques introduced in [64] is now possible to reach an higher precision and hence it 
will be possible to perform stringent tests, also on the spectrum of excited states. 

The second part of this work was dedicated to the study of scahng and universal- 
ity properties of gauge actions alternative to the plaquette formulation; in particular 
we concentrated our attention on the Iwasaki and DBW2 actions, which contain also 
rectangular (1x2) loops and have been formulated by using Renormalization Group 
arguments. 

Our main purpose was to test the universality and the scaling behavior for some 
relevant physical observable: we chose in particular the deconfinement temperature 
Tc and the mass of the glueballs 0^"*" and 2"^"''. The first issue that we faced was 
the setting of a suitable scale to build renormalized dimensionless quantities. As al- 
ready pointed out in previous studies it is important to use a quantity which can be 
measured with good precision and which is not subjected to uncontrolled systematic 
errors; the intermediate length tq had been introduced with this purpose and turned 
out to be a good choice. 

We performed new numerical simulations at several lattice spacings to compute ro/a, 
which was up to now not available on the literature for the Iwasaki and DBW2 ac- 
tions. In the analysis of our results we remarked that for these actions the physical 
positivity is violated [39]: in particular, negative contributions are expected in Eu- 
clidean correlation functions for small time. By investigating the location of the 
poles in the propagator at tree level, we evaluated for different actions the value of 
tjnin such that for t 3> tmin thesc unphysical contributions are expected to disappear. 
We discussed the implications for the extraction of physical observables, in partic- 
ular for the applicability of the variational method, which is founded per definition 
on physical positivity. 

Interesting features of these actions could also be learnt by studying the force at 
tree level: we noticed that RG actions are "overcorrected" at tree level and that 
beside the leading lattice artifacts of order a^, also higher powers are expected to be 
important. 

With our determinations of ro/a we then formed the quantity T^ro and we compared 
the results at different lattice spacings with the ones obtained with the plaquette 
Wilson action. For the Iwasaki action we observed a reduction of the lattice arti- 
facts, while for the DBW2 action the situation is less clear since the results depend 
quite strongly on which definition of the force at finite lattice spacing is adopted. 
The important result is that universality is confirmed for this particular observable; 
in previous calculations the quantity T^/ \fa was used and a discrepancy in the con- 
tinuum extrapolations between Iwasaki and Wilson actions was observed [38]. Our 
result exclude a violation of the universality and seem to confirm that this disagree- 
ment was generated by the problems related with the extraction of the string tension 
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a and hence by the choice of the quantity used to set the scale. 
With our results we were able to give a parametrization formula which describes ro/a 
in the range 2.1551 < /3 < 2.7124 for the Iwasaki action and 0.75696 < /3 < 1.04 for 
the DBW2 action, which can be useful for future applications. 

We checked the scaling of RG actions also on the running coupling extracted from 
the force by comparing the results at our smallest lattice spacing with the contin- 
uum results obtained in the first part of the work. As expected, at sufficiently large 
distances, we found no significative deviation, while at small r/ro we observed a 
difference of 2% for the Iwasaki action and 4 — 10% for the DBW2 action, depending 
also in this case on the definition of the force at finite lattice spacings. 
Finally we computed the glueball masses for the 0^+ and 2++ states. We started by 
selecting a basis of "good" operators from Wilson loops up to length eight, using the 
signal/noise ratio as criterion for the choice. 

For the evaluation of the glueball masses the violation of physical positivity is even 
more problematic that for the potential, since one can not avoid to use variational 
techniques. We analyzed our data by truncating the correlation matrices to suitable 
subspaces where the presence of unphysical states is supposed to be reduced. We 
build the quantities romo++ and rom2++; our final results are however affected by 
large errors; additional studies and numerical techniques are necessary to deduce 
clear conclusions on the lattice artefacts for these observables. However, by consid- 
ering our determinations as upper limit for the glueball masses at different lattice 
spacings, one could conclude that the RG actions are not able to significantly reduce 
the discretization errors for mo++ with respect to the plaquette action. 
Further investigations on the properties of RG actions could be very useful to clarify 
this question more definitively as well as further open questions, also in view of the 
future unquenched simulations. 
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Appendix A 
Numerical results 



In this appendix the numerical resuhs for potential and force at finite lattice spacing 
are collected, for different gauge actions. 

A.l Wilson action 

The tables A. 2 and A. 4 report the data for potential and force at finite lattice 
spacing (chapter 2) for the Wilson action. For the couphngs 5.7 < /3 < 6.4 we took 
the data from Guagnelli et. al. [73], while in the range 6.57 < (3 < 6.92 the results 
were obtained from our new numerical simulations, rj and di define the tree-level 
improved observables and are expressed in eq. (2.52), eq. (2.84). 

A. 2 RG actions 

The tables A. 5 and A. 6 report the data for the potential and the force evaluated 
at finite lattice spacing with Iwasaki and DBW2 actions, which has been used in 
chapter 3. Also in this case ri defines the tree-level improved force. For the potential 
only the naive definition is reported. 
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Appendix A Numerical results 



Table A.l: Force and potential from the data of Guagnelli et al. [73] 



(5 ri/a F{ri) di/a 



a Vi{di) 



5.7 



5.95 



6.07 



6.2 



2.277 0.2160(14) 

3.312 0.1847(25) 

4.359 0.1730(43) 

2.277 0.11212(20) 

3.312 0.08319(27) 

4.359 0.07164(36) 

5.393 0.06613(48) 

6.414 0.06296(56) 

7.428 0.0606(12) 

2.277 0.09211(11) 

3.312 0.06427(13) 

4.359 0.05301(26) 

5.393 0.04771(28) 

6.414 0.04468(27) 

7.428 0.04262(38) 

8.438 0.04215(45) 

9.445 0.04087(73) 

2.277 0.07804(11) 

3.312 0.05135(14) 

4.359 0.04054(13) 

5.393 0.03511(20) 

6.414 0.03238(20) 

7.428 0.03018(25) 

8.438 0.02884(26) 

9.445 0.02813(27) 

10.451 0.02766(30) 

11.155 0.02752f3n 



1.855 

2.889 
3.922 
4.942 

1.855 
2.889 
3.922 
4.942 

5.954 
6.962 
7.967 

1.855 
2.889 
3.922 
4.942 
5.954 
6.962 
7.967 
8.971 
9.974 

I. 855 

2.889 
3.922 
4.942 
5.954 
6.962 
7.967 
8.971 
9.974 
10.977 

II. 979 



0.80463(94) 

1.0206(22) 
1.2053(44) 
1.3784(77) 

0.61794(16) 
0.73006(33) 
0.81325(55) 
0.88489(84) 
0.9510(12) 
1.0140(17) 
1.0728(30) 

0.571729(97) 
0.66384(19) 
0.72811(30) 
0.78116(55) 
0.82887(72) 
0.87355(89) 
0.9162(11) 
0.9583(13) 
0.9992(17) 

0.533457(82) 

0.61145(17) 
0.66279(28) 
0.70333(38) 
0.73844(54) 
0.77082(69) 
0.80100(85) 
0.8298(10) 
0.8580(12) 
0.8856(14) 
0.9131 f 16) 
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Table A. 2: Force and potential from the data of Guagnelli et al. [73] 



(5 Til a 


F(ri) 


d\/a 




6.4 




1.855 


0.488379(48) 


2.277 


0.064318(51) 


2.889 


0.552697(81) 


3.312 


0.039580(70) 


3.922 


0.59228(13) 


4.359 


0.029360(79) 


4.942 


0.62164(19) 


5.393 


0.024367(84) 


5.954 


0.64600(24) 


6.414 


0.02145(12) 


6.962 


0.66744(40) 


7.428 


0.01939(16) 


7.967 


0.68683(50) 


8.438 


0.01819(18) 


8.971 


0.70502(66) 


9.445 


0.01757(17) 


9.974 


0.72258(79) 


10.451 


0.01677(17) 


10.977 


0.73936(84) 


11.455 


0.01651(15) 


11.979 


0.75586(94) 


12.459 


0.01609(17) 


12.980 


0.7720(11) 


13.462 


0.01616(29) 


13.982 


0.7881(11) 


14.465 


0.01564(18) 


14.983 


0.8038(12) 


15.467 


0.01513(39) 


15.984 


0.8189(14) 
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Table A. 3: Force and potential in the short distance region. 





n/a 


F(ri) 


di/a 


a Vi{di) 


6.57 






1.855 


0.457898(71) 




2.277 


0.056623(65) 


2.889 


0.51452(12) 




3.312 


0.033391(97) 


3.922 


0.54795(21) 




4.359 


0.023752(94) 


4.942 


0.57170(28) 




5.393 


0.01904(12) 


5.954 


0.59074(35) 




6.414 


0.01629(13) 


6.962 


0.60703(41) 








±.000 


4391 R(A1 ) 




9 977 


u.uozouo^oo J 








O.OIZ 


U.UOUl ( ^^^UU J 


'? Q99 






4.359 


021054('75) 


4.942 


54267fl8) 




5.393 


0.016439(72) 


5.954 


0.55911(23) 




6.414 


0.013728(82) 


6.962 


0.57284(28) 




7.428 


0.012036(96) 


7.967 


0.58487(35) 




8.438 


0.010869(82) 


8.971 


0.59574(41) 




9.445 


0.010123(97) 


9.974 


0.60587(46) 


6.81 






1.855 


0.422617(22) 




2.277 


0.048833(26) 


2.889 


0.471411(47) 




3.312 


0.027650(32) 


3.922 


0.499062(66) 




4.359 


0.018860(31) 


4.942 


0.517921(90) 




5.393 


0.014471(32) 


5.954 


0.53239(11) 




6.414 


0.011870(40) 


6.962 


0.54426(13) 




7.428 


0.010163(54) 


7.967 


0.55437(18) 




8.438 


0.009072(56) 


8.971 


0.56344(21) 




9.445 


0.008267(52) 


9.974 


0.57171(24) 




10.451 


0.007701(58) 


10.977 


0.57941(27) 




11.455 


0.007232(62) 


11.979 


0.58664(30) 
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Table A. 4: Force and potential in the short distance region. 





n/a 


F{ri) 


di/ a 


aFi(di) 


6.92 






1.855 


0.408642(19) 




2.277 


0.046081(19) 


2.889 


0.454723(34) 




3.312 


0.025696(24) 


3.922 


0.480418(47) 




4.359 


0.017266(29) 


4.942 


0.497662(75) 




5.393 


0.012969(33) 


5.954 


0.510631(88) 




6.414 


0.010412(39) 


6.962 


0.52104(11) 




7.428 


0.008855(37) 


7.967 


0.52990(13) 




8.438 


0.007755(43) 


8.971 


0.53765(16) 




9.445 


0.006974(51) 


9.974 


0.54463(19) 




10.451 


0.006402(53) 


10.977 


0.55103(23) 




11.455 


0.006022(53) 


11.979 


0.55705(26) 



Table A. 5: The potential and the force in lattice units for the Iwasaki action. 



/3 


ri/a 


a^F{ri) 


r/a 


aV{r) 


2.1551 


1.4858 


0.4059(21) 


2 


0.9519(20) 




2.4626 


0.2965(75) 


3 


1.2487(61) 




3.5663 


0.2585(98) 


4 


1.515(12) 




4.6059 


0.249(29) 


5 


1.783(12) 


2.2879 


1.4858 


0.6055(29) 


2 


0.7730(12 




2.4626 


0.19796(76) 


3 


0.9707(33) 




3.5663 


0.1689(70) 


4 


1.1418(41) 




4.6059 


0.1622(64) 


5 


1.3047(95) 




5.6082 


0.155(10) 


6 


1.460(16) 




6.6036 


0.148(16) 


7 


1.608(33) 


2.5208 


1.4858 


0.21301(25) 


2 


0.60342(29) 




2.4626 


0.11734(39) 


3 


0.72099(46) 




3.5663 


0.09011(70) 


4 


0.8107(14) 




4.6059 


0.0803(10) 


5 


0.8902(20) 




5.6082 


0.07575(61) 


6 


0.9658(36) 




6.6036 


0.0744(32) 


7 


1.0378(37) 




7.5977 


0.0711(15) 


8 


1.1098(74) 


2.7124 


1.4858 


0.17512(12) 


2 


0.52165(14) 




2.4626 


0.08691(14) 


3 


0.60856(20) 




3.5663 


0.06138(53) 


4 


0.66994(56) 




4.6059 


0.05194(81) 


5 


0.7219(13) 




5.6082 


0.04693(85) 


6 


0.7688(20) 




6.6036 


0.0451(23) 


7 


0.8139(11) 




7.5977 


0.0425(10) 


8 


0.8564(10) 




8.5915 


0.0420(15) 


9 


0.8984(23) 




9.5851 


0.0 11 If 19) 


10 


0.9398f25) 
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Table A. 6: The potential and the force in lattice units for the DBW2 action. 



p 


ri 


a r [Ti ) 


r /a 




U. ( obyb 




n ooo /'1 o\ 
O.oooo(^io ) 


o 

z 


n nnn /i /' 1 '7\ 
0.9004(^1 / j 




Z. ( /yd 


n oToo/'c^nN 
O.z/zo(oyj 


o 
o 


i.i /zb(33j 




Q '7/1 Q C 


n o/i 1 ^"10^ 


4 


1 AT T \ 
i.4iU(^ii J 




4. / C)04 


U.z4z4(^oo J 


e. 



i.bool^zU ) 




0. /OoO 


n oofi/'on'\ 


P. 

b 


i.oy4(^oi3 ) 


n oo /I o 


1 nooo 


n on^^n'y/^^ 1 \ 
U.zyby / (Oij 


o 
Z 


n ^70 ^7/ 1 1 \ 

0. / 3o j 




o TTno 
z.77yo 


n 1 ooT/i o\ 
0.1887(lzj 


o 
3 


n no / 1 1 \ 
0.yz44(llj 




O ^7/1 O 


n i^^io/'io\ 
U.ibio^iz j 


4 


1 no 1 
l.Uoo /(1 1 ) 




/I ^70^^ /I 

4. i oo4 


n T F;oo/'^^o^ 
U.iOoZ(^bZ j 


D 


1 ooon/'vn^ 

i.zooyi^ / U ) 




0. /ooo 


n 1 /I ^7^^ /n 1 \ 
U.i4/ bl^yij 


b 


l.oob(lb ) 




^70 /I o 


U.i4y4(bij 




1 txo^^/'on\ 
l.oob^zO j 


(J.yDOD 


1 nooo 

i.yzoo 


n onoin/'o'7\ 

0.z0810(z7j 


o 
z 


0.55by3(b4j 




o TTno 


0.113b(lbj 


o 
3 


0.6703(20) 




O. ( 4o0 


U.UoboU(^oD ) 


/I 
4 


n 7c;«n/'on'\ 
U. ( oVi)[A\J ) 




4.7364 


0.0783(13) 


5 


0.8340(29) 




5.7535 


0.0745(37) 


6 


0.9076(57) 




6.7843 


0.0700(57) 


7 


0.977(10) 




7.8120 


0.0640(47) 


8 


1.038(10) 


1.04 


1.9233 


0.18182(26) 


2 


0.49710(24) 




2.7793 


0.09392(54) 


3 


0.59094(58) 




3.7435 


0.06847(84) 


4 


0.65914(76) 




4.7364 


0.0575(18) 


5 


0.7180(26) 




5.7535 


0.0544(12) 


6 


0.7724(21) 




6.7843 


0.0527(22) 


7 


0.8254(42) 




7.8120 


0.0498(23) 


8 


0.8753(65) 




8.8296 


0.0481(42) 


9 


0.9230(51) 




9.8369 


0.0485(44) 


10 


0.9713(56) 




10.8365 


0.0464(34) 


11 


1.0175(95) 




11.8310 


0.048(13) 


12 


1.067(20) 



Appendix B 



Transfer matrix and Hamiltonian 
formalism 



The Hamiltonian formalism constitutes a very useful framework in order to extract 
physical observables from the exponential decay of Euclidean correlation functions. 
In this appendix the transfer matrix formulation will be reviewed, with a particular 
attention to the improved actions, where the violation of physical positivity at cutoff- 
energies can be a problem for the extraction of physical observables. 
In the formulation of eq. (1.8), lattice gauge theories are represented as a statistical 
mechanical system with partition function Z and thermal averages (C). The transfer 
matrix provides an equivalent representation of the model as a quantum mechanical 
system with a Hilbert space 7i, a Hamiltonian operator H and linear operators O 
corresponding to the observables O. 

We will now assume that the Euclidean space-time lattice A is finite, with T x L x 
L X L points and we take periodic boundary conditions for the gauge fields 



In this appendix we will set a = 1 for simplicity. 

B.l Wilson action 

We denote the set of links that are completely contained in the time slice xq = t hy 
Bt] those links which connect the time slices xq = t+1 and xq = t form the set Bt+i,t- 
Let 6 be a link connecting the lattice points x and {x + fi) and U{h) = U (x, //.); then 
the links belonging to Bf and Bt+i^t are denoted by 



U{x + aTu, n) 
U{x + aLi), /i) 



f/(,T,/i), z/ = 

U{x,iJ,), u > 0. 



(B.l) 
(B.2) 



Ut = U{h), beB, 

Ut+i,t = U{b), beB, 



't+l,t- 



(B.3) 
(B.4) 



The Wilson action can be written as 



S = J2L[Ut+i,Ut+i,t,Ut], 
t 



(B.5) 
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where 

L [t/t+i, Ut] = ^Li [[/i+i] + ^Li [Ut] + L2 [C/t+i, C/i] . (B.6) 

The first two terms in eq. (B.6) 



a; fc<i 



are the contribution of space-hke plaquettes on a time shce, while 

L2[t/,+i,t/,+i„t/,] = -^Y.Y.^''^lo\^)U=t (B.8) 

X k 

represents the contribution of time-hke plaquettes between two time slices. The 
plaquette W^^^{x) is defined in eq. (1.4). 

We then consider wave functions ^[C/*] depending on the link variables on a certain 
time slice; these describe states belonging to different charged sectors depending on 
how they transform under gauge transformations. In particular, we will consider 
gauge invariant wave functions, which correspond to the vacuum sector 

= (B.9) 

where ^7^ and Ut are gauge-equivalent. 

A convenient way to define a scalar product of two given wave functions ^ and $ is 
given by 

($,^)= J DptMUtYnUt], (B.IO) 

with 

D[Ut] = n dU{b). (B.ll) 

beBt 

The Hilbert space H corresponds to the space of all gauge invariant functions with 
finite norm, 

(*,*)< 00. (B.12) 
The transfer matrix T [131] is an integral operator acting on Ti via 

(T*)[C/t+i] = J D[Ut]K[Ut+i,UtMUt], (B.13) 
where K is an integral kernel 

K[Ut+i, Ut] = J D[Ut+i,t] exp{-L [Ut+i, Ut+i,t, Ut]}. (B.14) 
Prom these considerations follows that 

J D[U]e-^^^^ = tr {T^}, (B.15) 
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where "tr" indicates the trace in the Hilbert space. 

The Euchdean correlation functions {A{x)B{y)...) of local, gauge invariant fields 
A{x), B{y)... can be given in a quantum mechanical interpretation. In particular, if 
these fields are polynomials of the gauge field variables Ut at a fixed time we may 
associate an operator A{x) through 

{A{x)-^)[Ut]^ A{Q,x)-^[Ut]. (B.16) 

For = t, < < T we have 

{A{x)A{Q)) = I tr {T^-*i(f )T*i(0) } . (B.17) 

Similar formulas hold for higher correlation functions. The transfer matrix T turns 
out to be a bounded, self-adjoint positive operator with a completely discrete spec- 
trum (for finite L). 

All eigenvalues of T are then strictly positive and the maximal eigenvalue Aq is not 
degenerate. Thus one can define the Hamiltonian 

T 

H=-ln— , (B.18) 

which is a self-adjoint operator in V,. For T — > oo the two-point function eq. (B.17) 
then assumes the form 

{A{x)A{{))) = (0|i(f)e-*^i(0)|0), (B.19) 

which is also called spectral representation and establishes the exponential decay 
of Euclidean two-point functions, which is the basis of the extraction of physical 
observables from lattice QCD. 



B.2 Improved actions 

We now consider actions containing not only the plaquette term, but also other 
classes of loops C included in a real gauge invariant and orientation independent 
function of the ordered product of link variables along C. These include both 
Liischer-Weisz and RG improved actions as particular cases. One can prove that 
for those actions it is still possible to construct a transfer matrix [39] . Restricting to 
actions of kind eq. (1.13), one can notice that in the temporal gauge, C/(x, 0) = 1, 
the associated field equations assume the form of a fourth-order difference equation 
in the time coordinate Xq. One has then to define a transition amplitude to go from 
a given gauge configuration on a pair of consecutive equal time hyperplanes to some 
other gauge fields on a later pair of hyperplanes. For the wave function ^ one then 
makes the ansatz 

* = C/t+i,i,C/t]. (B.20) 

A gauge invariant scalar product is now 



($, *) / D[Ut+i]D[Ut+i,t]D[UMUt+uUt+i,u Ut]*-^[Ut+i, Ut+i,t, U^]. (B.21) 
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The transfer matrix can be defined in analogy to the Wilson case as transition am- 
plitude between fields configurations defined on two consecutive pairs of hyperplanes 
{Bt,Bt+i) and (5t+i,5t+2), and one can argue that the properties eq. (B.15) and 
eq. (B.16) are satisfied. 

But, unlike the Wilson case, one finds that T is no longer hermitian; rather, one 
obtains 

T^" = eje, (B.22) 

where 

9^Ut+^,Ut+^,t,Ut]^^Ut,u^+%Ut+,], e^^e, e^^i. (B.23) 

The spectrum (t(T) of T can be defined as the set of complex numbers A, for which 
the operator A — T has no bounded inverse. Further considerations [39] lead to the 
conclusion that spectral values always occur in pairs of complex conjugated numbers 
and there is a positive (non-degenerate) eigenvalue Aq such that 

Ao>|A|, VAe(7(T), A^Ao. (B.24) 

The eigenfunction ^x^^ may be interpreted as ground state wave function. Due to 
the violation of physical positivity, it is no longer possible to define an Hamiltonian 
operator and hence to write down the spectral representation eq. (B.19). For T ^ oo 
the two-point function eq. (B.17) takes the form 

{A{x)A{0)) = Ao*(0|M(f)T*i(0)|0), (B.25) 
where |0) = I^'ao); a-iid in the limit t ^ oo one ends up with an asymptotic expansion 

(A(x)A(O)) ~*_oo Ao* Yl 

\ea{T)\\o 

where px{t) is a polynomial of t with degree strictly less than a certain nx: usually 
nx — 1 as in the hermitian case, but one can find examples where nx > 1 for 
some A. The eq. (B.26) establishes that connected two-point functions always decay 
exponentially at large times, but if the leading spectral value is not real and positive, 
the exponential factor is multiplied by an oscillating amplitude and contributions in 
the spectral decomposition with negative weight appear, as remnant of the positivity 
violation. This fact has been observed in the evaluation of the effective potential 
and effective masses for improved actions 3, 4 and can be a serious problem for the 
applicability of the variational method. 

Nevertheless, one expects that physical positivity is recovered in the continuum limit; 
restricting ourself to the subspace Hphys of the states with small energy with respect 
to the cutoff (A/Ao ~ 1), then the non-hcrmiticity of the transfer matrix can be 
eliminated by choosing a new scalar product relative to which T becomes hermitian 
and physical positivity is then completely restored. For example 

{^,^)new^i^,0'^)old (B.27) 

is not well defined on the full Hilbert space, but one expects that there exists a 
< e < 1 such that independently of the cutoff the following properties hold: 
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(i) all spectral values A G (t(T) with |A| > eAo are real and positive; 

(ii) ^)nev. > 0, V* G Hphys, ^ 7^ 0. 

The existence of such e has not been rigorously proved, but is supported by pertur- 
bativc calculations, which can be used to estimate a value of e by determining the 
location of unphysical poles in the propagator, as shown in the next section. 

B.2.1 Unphysical poles in the propagator 

In [47] the propagator D/^i, is evaluated to lowest order in qq with covariant gauge 
fixing; it is defined by 

{A'^{x)Ai{y)) = 5'^ ! e'''^^-y^e'^'--'^y^D^,{k), (B.28) 

Jk 

where Aj^i^x) is related to the link variables by 
and the momenta integration corresponds to 



k 



We arc now referring to the case L = oo, so that the momenta take continuous 

values in the Brillouin zone. 

Eq. (B.28) can be rewritten in the form 



D,r{k) = {k 



2\-2 



(B.29) 



where a is the gauge parameter, A;^ = 2sin(/c^/2) and A^^, is independent of a and 
satisfies the following properties: 

1. A,^^ = 0, V//; 
9 A — A ■ 

3. A^^i^k) = y4^i/(— /c); 

4. A^,{0) = 1, ^ I/. 
A^i, has the general form 



A,. = (B.30) 



where 



<lij.v + k'^J^'l(lfj.v{<liJ.p<lvT + %t(1up)-, (B.31) 

H Vi^p, //>I/,P>T,{/),t}U{;U,I/}=0 
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and the expUcit form of f{k) can be found in [47] and will not reported here. The 
denominator of eq. (B.29) is then given by 

A = D{kY = (B.32) 

\ M J'T^M M>'^,/O>-r,{p,-r}U{/i,i/}=0 

where q^v satisfies the properties (l)-(4) and in our specific form of the action takes 
the form 

q^, = {l-5^,)[l-c^{kl + kl)], (B.33) 

Ci being the coefficient of the 1x2 planar loops term in the action eq. (1.34). Then 
eq. (B.32) becomes 



A= p-cj^fc; 



P-c,(^(P)^ + 5^^jj+ (B.34) 



In order to search for the poles of the propagator, one substitutes 

k = {ki,k2,k3,iw) 

and looks for solutions of the equation A = scanning the whole Brillouin zone. 
For the plaquette action ci — one finds the usual result 



w = ±\/kl + A;| + fc| 

while for general Ci one expects complex conjugated solutions 

w = ^t{w) ± iJm(w) = w{ki, k2, ka). 

Numerical investigations showed that the condition '3m{w) = (for all solutions 
at a given momentum) defines the equation f{ki, k2, k^) = which represents a 
compact 3-dimensional object. In Fig. B.l the 2-dimensional intersection with the 
plane ^3 = is plotted for several values of ci corresponding to different actions. 
For our numerical studies, we discretized the Brillouin zone in finite intervals that 
could be made arbitrarily small. 

In the region enclosed by the curve one has 3m{w) = 0, while outside 3m{w) 7^ 0. 
The Fig. B.2 represents the distribution of the poles; by scanning the Brillouin zone 
and varying the intervals of the momenta, one reaches the conclusion that there 
exists a maximal value of ^Rc{w), which is related to Ci by 

2.063 Symanzik, tree level 
Wmnr = 2arcsinh ( — . 1 = i 1.162 Iwasaki 
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Figure B.2: Distribution of the poles for different actions, scanning the Brillouin 
zone. 
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such that for w < Wmax the imaginary part 3m{w) vanishes. 

The value of Wmax yields an estimate of e defined above through the simple relation 

{0.127 Symanzik, tree level 
0.313 Iwasaki 
0.555 DBW2 

Evaluating observables from exponential decays of correlation functions, one expects 
the presence of unphysical states unless t ^ tmim where tmin in lattice units is given 

by 

^ ( 0.484 Symanzik, tree level 

tmin — = { 0.860 Iwasaki (B.35) 

Wmax y ^ 702 DBW2 

This fact has to be taken into account in the extraction of physical observables from 
correlation functions evaluated with RG actions, and in general with every action 
which contains not only plaquette terms. In particular this restricts the applicability 
of the variational method. 

One has to remember that our estimation of the unphysical poles has been performed 
in perturbation theory and one can not exclude a priori that at large coupling qq the 
numerical value of e can be larger than what we found. 



Appendix C 

Evaluation of the 3-dimensional 
lattice propagator in coordinate 
space 

An efficient metfiod to calculate the lattice propagator in coordinate space is pro- 
posed in [132]. It is based on a recursion relation which allows to express the propa- 
gator as linear function of its values near the origin. In that paper, the 4-dimensional 
case is discussed; it is straitghforward to apply the same method in 3 dimensions. We 
use lattice units, a = 1, in tliis appendix and restrict the attention to the plaquette 
Wilson action. 

C.l Lattice 3-dimensional propagator. 

The 3-dimensional scalar propagator in the coordinate space satisfies the Laplace 
equation 

^G{x) I Q otherwise. i^-^) 
where the lattice Laplacian is given by 

and X = {xi, X2, x^). 

The forward and backward lattice derivatives are defined through 

Vjf{x)=f{x + j)-f{x), V*f{x)=f{x)-f{x-j), (C.3) 
where j denotes the unit vector in direction j. 

In the continuum limit, which here coincides with the limit \x\ oo, we expect that 
the propagator converges to 47r|a;|, which is the Green function of the continuum 
Laplacian on R^. 



108 Appendix C 3- dimensional lattice propagator in coordinate space 



C.2 Recursion relation 

The first step to obtain a recursion relation is the observation that 

{V* + V,)G{x)=x,H{x), (C.4) 

where 

is independent of the direction j and 

3 

P = ^ A?/, kj = 2 sin(%/2). (C.6) 
Summing over k and using eq. (C.l) one obtains 

H{x) = -J2[G{x)-G{x-j)]. (C.7) 



with p = Yl^.-^ 



j=i -J 



This expression can be used to eliminate H(x) in eq. (C.4), giving 



X 



G{x + j) = G{x - j) + 2^ ) - G{x - i)] (C.8) 

p tr 

which is valid only for the points with p 7^ 0. 

Because of isotropy, we can restrict our attention to the points x with xi > X2 > 
> 0. In this region eq. (C.8) is a recursion relation which can be used to express 
G{x) as a linear combination of the values of the Green function at the corners of 
the unit cube, 



^(0,0,0), ^(1,0,0), ^(1,1,0), ^(l,!,!). (C.9) 
This four initial values are not independent; from eq. (C.l) at x = 

G(0,0,0)-G(1,0,0) = ^ (C.IO) 

o 

follows directly. Another relation can be deduced in the following way: first one can 
observe that eq. (C.8) becomes one-dimensional along the lattice axes. Defining 

^i(n) = G'(n,0,0), ^2(n) = G(n, 1, 0), g^{n) ^ G{n,l,l), 

and using the lattice symmetries, one finds 

^i(n+l) = 6^i(n)-4^2(n)-^i(n-l), (C.ll) 

2?7/ Tt — 1 

^2(^+1) = ^-[3^2(n)-^i(n)-y3(n)]-^-^72(n-l), (C.12) 



C.3 Numerical computation of G{x) 
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2n n — 2 

gsin + l) = [3g,{n)-2g2{n)]-^-gs(n-l). (C.13) 
n + 2 n + 2 

Next one notices that 

I{n) = (n - l)gi{n) + 2ng2{n) + {n + l)g3{n) - ngi{n - 1)- (C.14) 

2(n-l)(72(n-l)-(n-2)^3(n-l) 
is an invariant of the recursion relation, that is 

7(n + l) =7(n) 

for n> 1. 

The value of I can be worked out in the limit n — > oo, where 

gj{n)^^ + 0{l/n^), for j = 1,2,3, 
yielding 7 = 0. Setting n — 1, eq. (C.14) gives 

3^(1, 1, 0) + 2G(1, 1, 1) - ^(0, 0, 0) = 0. (C.15) 

Thus from the four initial values two can be eliminated using eq. (C.IO) and eq. (C.15) 
and the propagator can be written in the form 

G{x) = ri{x)G{0, 0, 0) + r2(f)G(l, 1, 0) + rs{x), (C.16) 

where the coefficients ri, r2, are rational numbers that can be computed recursively 
from eq. (C.8). 



C.3 Numerical computation of G{x] 

Prom eq. (C.16) one obtains the initial conditions for ri,r2,rs : 



ri(0,0,0) = l, ri(l,l,0) = (C.17) 
r2(0,0,0) = 0, r2(l,l,0) = l (C.18) 
r3(0,0,0) = 0, r3(l,l,0)=0. (C.19) 

Eq. (C.IO), eq. (C.15) give 

ri(l,0,0) = l, ri(l,l,l) = | (C.20) 
r2(l,0,0) = 0, r2(l,l,l) = -| (C.21) 

r3(l,0,0) = -^, r3(l,l,l) = 0. (C.22) 

The other coefficients ri{x) , r2{x) , rs{x) can be calculated exactly using for example 
an algebraic programming language. We are then left with the task to compute 
numerically G{0, 0, 0) and G{1, 1, 0) with high accuracy. For this purpose we consider 
the lattice points x — (n, 0, 0) and y = (n, 1,0), where n is an adjustable integer. 
The numerical investigations show that the associated coefficients rk{x) and rfe(y) are 
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Table C.l: Numerical values of the coefficients rk{x), rk{y) for n = 1, .., 7. 



n 


ri{x) 


ri{y) 


r2{x) 


r2{y) 






1 


1 








1 


-0.17 





2 


5 


-1.5 


-4 


4.5 


-1 


0.17 


3 


35 


-14 


-42 


28.3 


-6.5 


2 


4 


261 


-123.75 


-361.3 


219.75 


-46 


19.17 


5 


2026 


-1091.2 


-3005 


1844.413 


-346.17 


174 


6 


16259.8 


-9663.5 


-25046.32 


16000.9 


-2727 


1559.17 


7 


134186.8 


-86122.9 


-211276.52 


141288.85 


-22252.5 


13968 



rapidly growing, roughly proportional to 10" ^ (see Table C.l). Solving eq. (C.16) 
for G{x) and G{y), we obtain 

G(,0.0)^i{o(.,-G(.,^.Ii^-.3(.,}. (C.3) 



with 



C = r.fe)-r2MrM, (C.24) 

ri{x) 



and 



G(0, 0, 0) = -^G(l, 1, 0) - (C.25) 

The coefficients multiplying G{x) and G{y) are of order lO-^""^) and can be ne- 
glected, obtaining 

This method is exponentially convergent and the error can be estimated well since 

the neglected terms are known. 

By applying this technique one obtains 

^(0,0,0) = 0.2527310098586630030260020266135701299..., 
^(1,0,0) = 0.0551914336877373170165449460300639378... 
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